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DEVELOPMENT AND APPLICATION OF A METHOD FOR 


PREDICTING ROTOR FREE WAKE POSITIONS AND 
RESULTING ROTOR BLADE AIR LOADS 
VOLUME I - MODEL AND RESULTS* 

By S. Gene Sadler 

Rochester Applied Science Associates, Inc. 


SUMMARY 


Rotor wake geometries were predicted by a process similar to 
the start-up of a rotor in a free stream. An array of discrete 
trailing and shed vortices is generated with vortex strengths 
corresponding to stepwise radial and azimuthal blade circulations. 
The array of shed and trailing vortices is limited to an arbitrary 
number of azimuthal steps behind each blade. The remainder of the 
wake model of each blade is an arbitrary number of trailed vortices. 
Vortex element end points were allowed to be transported by the 
resultant velocity of the free stream and vortex-induced veloci- 
ties. Wake geometry, wake flow, and wake-induced velocity influ- 
ence coefficients were generated by this program for use in the 
blade loads portion of the calculations. 

Blade loads computations included the effects of nonuniform 
inflow due to a free wake, nonlinear airfoil characteristics, and 
response of flexible blades to the applied loads. The resulting 
nonlinear equations were solved by an iterative process to deter- 
mine the distribution of blade shears, bending moments and twisting 
moments . 

Computed wake flows and blade loads were compared with experi- 
mentally measured data. Predicted blade loads, response, and 
shears and moments were obtained for a model rotor system having 
two independent rotors. The effects of advance ratio, vertical 
separation of rotors, different blade radius ratios and different 
azimuthal spacing of the blades of one rotor with respect to the 
other were investigated. 

Wake flow calculations indicate that the wake geometry model 
provides a satisfactory simulation of wake-induced flow near the 
rotor disk. Free wake geometry calculations indicate increasing 
distortion with decreasing advance ratio, intertwining of pairs 
of tip vortices for rotor systems with nonuniform azimuthal blade 
spacing, and more vortex-blade interactions than would be expected 
from a rigid wake model. Blade airload variations are generally 
consistent with wake geometry calculations, but do not show ef- 
fects of vortex-blade interactions as might be displayed by a 
more detailed model. 

ic ’ * 1 
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INTRODUCTION 


The prediction of rotor-induced aerodynamic flows has recently 
undergone extensive refinement, primarily due to the need for im- 
proved predictions and the availability of high-speed large-capacity 
digital computers. The need for increased performance, lower noise 
levels, and other considerations have stimulated the development of 
improved analytical tools for use in helicopter design. The capa- 
bilities of modern computers together with the need for refined 
analyses has led to more realistically detailed analytical models 
for calculation of approximate nonuniform wake-induced inflow, 
blade natural frequencies and mode shapes, and blade loads. The 
improved approximations for the wake model have provided a fair 
quantitative measure of the nonuniform wake-induced inflow, and 
due to this , an improved prediction of the higher-harmonic air- 
loads acting on the blades. 

In recent times, the development of a model for nonuniform in- 
flow determination was carried out independently, and at about the 
same time by Piziali and DuWaldt (ref. 1) and by Miller (ref. 2) . 

The many applications of this theory, such as those of references 
3 and 4 , have made quite clear the importance of wake structure 
and geometry in determining blade dynamics. In particular, the 
outboard portion of the wake of each blade, which rapidly rolls 
up into a concentrated vortex (usually termed a tip vortex) has a 
significant influence on higher-harmonic loading (refs. 3 and 5), 
transient response (ref. 6) , rotor noise (ref. 7) , and overall 
performance (ref. 8) . 

Improved high-harmonic airloads alone are of limited useful- 
ness in response and blade stress analysis unless compatible blade 
models are used in blade response calculations. In order that the 
response to higher-harmonic airloads be adequately computed by the 
use of normal modes of free vibration as generalized coordinates, 
modes with frequencies corresponding to the larger range of air- 
load frequencies must be used. Current programs allow determina- 
tion of these quantities for fully coupled modes for rotating 
blades, so compatible models for blade loads and blade response 
are available. 

The determination of wake- induced flow at the rotor blades 
has been the most difficult part of the model to develop satis- 
factorily. This is primarily a result of the sensitivity of wake- 
induced flow to the magnitude of the relatively small distances 
between the tip vortex of one blade and the subsequent blade, 
especially on the advancing side. Similarly, because of the 
interactions of the various wake elements, the location of the 
inboard wake structure also has a significant effect on blade 
loading. Thus a reasonably accurate definition of the wake 
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geometry is needed in order that wake effects can be properly 
evaluated. 

A rigorous determination of the structure and motion of the 
wake of a rotor is an extremely arduous task. The flow is unsteady, 
three-dimensional, and allows no convenient basis for lineariza- 
tion. As a result, the wake has been assumed to be a fixed skewed- 
helical pattern for forward flight conditions, and wake distortions 
have generally been neglected. Studies reported in refs. 9 and 10 
indicated the possibility of predicting wake distortion effects. 

In reference 9, the flow due to a rotor hovering near a ground 
plane was approximated by the axisymmetric flow due to a succession 
of vortex rings released from the rotor plane. A direct time inte- 
gration of the flow was carried out, each vortex being free to con- 
vect at the local induced velocity. A fairly good approximation 
to the physical flow being modeled was obtained. There was clear 
evidence of wake contraction and the computed and measured flow 
exterior to the wake envelope agreed well. Near the rotor plane 
the interior flow was poorly predicted; it is believed that this 
resulted because the inboard wake structure was not represented. 

In reference 10, the flow due to a rotor in forward flight was 
calculated. The wake of each blade was represented by a single 
concentrated tip vortex having a finite core of rotational fluid. 

The wake location for each rotor azimuth was determined by ini- 
tially positioning the wake at locations approximating a skewed 
helix, allowing the wake to become free, and then performing con- 
tinuous time integrations until a periodic flow was established. 

The calculation of the velocity at which each wake element was 
convected at a given instant included a term depending on curva- 
ture and core size of the vortex to account for local self-induced 
effects. Good agreement with flow measurements some distance from 
the rotor plane was obtained. Again, however, it is believed that 
omission of the inboard wake structure caused large errors in pre- 
dicting the flow in or near the rotor plane. 

Some attempts have been made to account for wake distortion 
in blade load analysis. In the calculation of VTOL propeller per- 
formance in reference 11, determination of wake distortion was 
made a part of the overall calculation. In this wake, considera- 
tion was limited to the hover case which made the flow steady rela- 
tive to the blade. Reference 12 reports the results of the analysis 
of hovering performance of rotors as calculated using experimental 
results to position the wake. In a similar approach, the theoret- 
ical results of reference 9 were used in reference 8 to make an 
estimate of wake contraction for VTOL propellers in yawed flight. 
Thus, both theoretical and experimental wake distortion effects 
have been used in attempts to predict more realistic induced 
velocities. Although the procedures used in these studies are 
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limited to rather special cases, they do indicate the importance 
of accounting for wake distortions in the analysis of rotors and 
propellers . 

Calculation of airloads and blade response through higher 
harmonics is often vital to the analysis of blade shears and mo- 
ments, blade life and other design considerations. With the ad- 
vent of new rotor configurations, such as rotors of different 
physical properties or nonuniform azimuthal blade positioning, a 
program capable of analysis of such systems became desirable. The 
present program has the capability of predicting wake geometry 
and wake flow, including the induced flow in the rotor plane for 
use in blade loads calculations, for a system of one or two heli- 
copter rotors. If two rotors are used for the helicopter model, 
they may be arbitrarily located with respect to one another, have 
blades of different properties, rotate in the same or different 
directions, and the wakes from all blades are allowed to interact 
and to affect the induced velocity distribution on each of the 
rotor blades. The blade loads program uses input from both an 
independent blade frequency and mode shape prediction program and 
from the wake geometry program discussed herein to perform calcu- 
lations of blade loads and blade response by an iterative proce- 
dure. Upon termination of these calculations, blade shears and 
moments are computed. 

In addition to the calculation of distorted wake geometries, 
the effects of such geometries on wake-induced flow, and the result- 
ing blade airloads, shears and moments, the program may also be 
used to compute other quantities of interest. For example, rotor 
angle of attack distributions are now computed but not written out, 
as are approximate blade loads during the wake geometry calcula- 
tions. Similar predictions for models of helicopter rotor blade 
response in maneuvers, under gust loading, and similar applica- 
tions may be possible with modifications to the existing program. 

By using a starting rotor approach with forward time integration 
as the basic method for calculation of wake geometry and wake- 
induced velocity effects, it is expected that this approach and 
program may be extended to the analysis of non-steady-state flight 
conditions . 


SYMBOLS 


a . 


1 


core radius of i— vortex element 


A , A , A , mode shape quantities representing the torque, 

v w linear flap deflection, linear lead-lag deflec- 

<j> ' ip tion, angular torsional deflection, and angular 

lead-lag deflection, respectively. Torque dimen- 
sions in meter-newtons, linear deflection in 
meters, angular deflection in radians. 
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nondimens ionalized semichord lenqth 
chord length, meters 
section drag coefficient 

section lift coefficient 

dC £ 

slope of section lift coefficient, -5 

da 

section moment coefficient (about midchord) 

thrust coefficient, (aircraft weight) / (ttpQ 2 R 4 ) 

perpendicular distance from point of induced 
velocity to axis of vortex element, meters 
(see fig. 3) 

generalized force of the K— mode 

lumped aerodynamic forces acting at blade mass 
points, newtons 

aerodynamic forces per unit span acting tangen- 
tial and normal to the rotor plane, respectively 
newtons/meter 

mass moment of inertia of blade about flap hinge 
meter-newtons/second 2 

torsional mass moment of inertia about c.g. of 
th 

i — section of the lumped parameter blade model, 
meter-newtons/second 2 

flap spring stiffness, meter-newtons/radian 
lift per unit span, pounds/foot 
blade mass, kilograms 
th 

mass of i— section of lumped parameter blade 
model, kilograms 

aerodynamic flapping moment on the blade about 
the blade hinge, meter-newtons 

Mach number 

aerodynamic twisting moment about midchord per 
unit span, newtons 


5 


lumped aerodynamic twisting moment about elastic 
axis , meter-newtons 

number of azimuthal steps per revolution 

number of radial load points per blade 

number of radial load points on all blades 
number of azimuthal steps in the wake 

vortex induced velocity, meters/second 

vortex self-induced velocity, meters/second 


forcing functions corresponding to linear flap, 
linear lead-lag, angular twist and angular lead- 
lag motions, used in computing generalized 
forces 

radial coordinate, meters or dimensionless 
radial distance to blade chord, meters 
radial distance, meters 
position vector of vortex element 

position vector of point at which induced velocity 
is computed 

position vector, s -s 

-p - 


magnitude of s 


5.1 ' 


Ei 


radius of curvature 
time, seconds 
torque , meter-newtons 

velocity tangential to rotor plane, meters/second 
linear flapwise blade deflection, meters 
velocity normal to rotor plane, meters/second 
aircraft forward velocity, meters/second 

wake-induced velocity, meters/second or linear 
edgewise blade deflection, meters 



the chordwise average downwash, meters/second 

the chordwise variable downwash, meters/second 

uniform wake-induced inflow, meters/second 

wake-induced inflow due to wake elements having 
known circulation, meters/second 

spanwise distance to blade center from flap hinge, 
meters 

tilt of tip path plane due to blade flap, posi- 
tive forward, radians 

sum of collective and cyclic pitch and blade 
twist angles 

shaft tilt angle, positive aft, radians 

blade flapping angle, B = Bo + Bi sin ftt 
+ 3 2 cos fit, radians 

steady and first harmonic lateral and longitudinal 
blade flapping components , radians 

circulation distribution , meters/second 

circulation distribution corresponding to uniform 
downwash, meters/second 

circulation distribution corresponding to variable 
downwash , meter s/second 

vortex element circulation, meters 2 /second 

circulation, chordwise integral of Yq, meters 2 / 
second 

circulation, chordwise integral of Yi, meters 2 /second 

distance between flap hinge and shaft axis, meters 

chordwise separation of the elastic axis and cen- 
ter of mass, positive forward, meters 

th 

generalized coordinate for the K — mode 

angles used in vortex-induced velocity determina- 
tion, see fig. 3, radians 



e 

c 

©o 

©i 

y 

Z 

p 

a 


$ 

$ 

¥ 


ft 


cyclic pitch control input angle, radians 

rigid body pitch, radians 

blade twist, radians 

advance ratio, V^/ftR 

local blade spanwise slope, radians 

air mass density, kilograms /meter 3 

wake-induced velocity influence coefficient 

blade elastic torsional deflection, radians 

angle between chord and rotor plane, radians 

azimuthal coordinate, positive counterclockwise, 
zero aft, radians 

blade elastic edgewise bending angle, radians 

th 

natural frequency of the K— mode, radians/ 
second 

rotor speed, radians/second 


Subscripts and Mathematical Terminology 


i , l , n 

j , k 

K 

v 


w 


X 


indicate radial station of variable 
indicate azimuthal station of variable 
indicates vibration mode number 

indicates force component acting in direction of 
v-displacement 

indicates force component acting in direction of 
w-displacement 

indicates force acting parallel x-axis 


z indicates force acting parallel z-axis 

cj) indicates moment component acting in direction 

of <j)-displ acement 
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jp indicates moment component acting in direction of 

^-displacement 

A indicates an increment, i.e. At is an increment 

in time 

. indicates time derivative, i.e. v = g— 


ANALYSIS 

The overall model and program arrangement is shown in Figure 
1. The model and problem formulation may be conveniently thought 
of in two sets, one dealing with the wake model and associated cal- 
culations, the other dealing with determination of the aerodynamic 
blade loads and response. The combination of these analyses and 
corresponding computer programs are directed toward the prediction 
of blade loads and associated dynamic response for helicopters in 
steady forward flight in which the effects of a free wake and flexi 
ble blades are included. One or two rotors may be modeled, and if 
two rotors are used, the blades on one rotor may differ physically 
from the blades on the other rotor. Blade-wake interactions are 
allowed, but no mass or elastic couplings between blades on a rotor 
or between blades on different rotors are allowed. While the num- 
ber of blades per rotor is arbitrary (subject to practical limits 
of the dimensioned variables in the computer program) , each rotor 
is assumed to have the same number of blades, and to have the same 
rotor speed. 


Wake Geometry Model and Formulation 

The wake geometry is calculated by carrying out a process 
similar to the startup of a rotor in a free stream. The blades are 
located at specified azimuthal and flapping positions, without any 
wake vortices. The blades then rotate through an azimuthal incre- 
ment, A\p , and shed and trail vortex elements of unknown strength, 
but with knov/n positions. The strengths of the vortices that are 
shed immediately behind the blade are then determined, and include 
the effect of their own self-induced velocities. (If it is de- 
sired, an estimate of the blade loads that result can then be 
determined, without the effects of blade flexibility being includ- 
ed.) All vortex element end points not attached to the blade are 
then allowed to translate as the blade is stepped forward for a 
time At with velocities as determined by the free stream and in- 
duced velocities. Here, At = Aij;/ft, where ft is the rotational speed 
This completes a typical first step in the wake geometry calcula- 
tion. 
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II 


Subsequent steps are similar. The blades are rotated through 
an angle and to new (specified) flapping angles. The vortices 

at and just behind the blades have unknown strength, while those 
in the wake have known strength. The strengths of the vortices 
attached to the blade -are then computed, induced velocities at 
vortex element end points are computed, and end points not at- 
tached to the blade are allowed to translate with the resultant 
velocity due to the free stream and induced velocities. (Again, 
if desired, approximate blade loads may be computed, once blade 
circulations are known.) In this manner arrays of discrete shed 
and trailing vortices are generated immediately behind the blades 
with strengths which correspond to approximate blade loads. These 
arrays have stepwise radial and azimuthal strength variations so 
that total circulation is conserved. The arrays of shed and trail- 
ing vortices which are generated immediately behind the blades are 
referred to as the full mesh wake. Comparisons of wake flows pre- 
dicted using this wake model for the entire wake with experimental 
measurements indicated that retention of shed elements with a 
coarse mesh resulted in poor induced velocity predictions , and 
that use of a fine mesh increased running time to an unacceptable 
level. Therefore, the full mesh wake was used to represent the 
wake immediately behind the blades, and a modified wake model was 
developed and implemented for use in the representation of the re- 
mainder of the wake, as shown in Figure 2 . 

The modified wake consists of trailing vortices only, so vor~ 
ticity is not conserved. The locations, circulations, core sizes, 
etc. of the trailing vortices in the modified wake are determined 
in terms of the full mesh wake's last sets of trailing vortex 
locations, circulations, etc. The wake- induced velocities, wake 
distortions, and other calculations are essentially the same for 
both the full mesh wake and for the modified wake portions of the 
wake model . 

The wake and blade vortices are assumed to have finite-sized 
cores of rotational fluid, and the core sizes at the blade are con- 
trolled by an input parameter so as to be adjustable to improve 
agreement between calculated and experimental or theoretical data 
as desired. All elements are assumed to be straight for purposes 
of calculations, except where induced velocity calculations are 
performed for an end point of a vortex element due to an element 
with one end at that point. At such a point the neighboring ele- 
ments are assumed to be arcs with curvature determined by three 
appropriate end points. 

The basic mathematical formulations of the wake geometry pro- 
gram are those related to determination of blade flapping motion, 
calculation of vortex- induced flow, and calculation of blade cir- 
culations including self-induced effects. 
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Determination of the approximate blade rigid-body flapping 
motions . - If blade flapping motions, collective and cyclic pitch 
angles, and flight conditions are known, these are input directly 
to the program. If flight conditions and collective and cyclic 
pitch angles are known but blade flapping motions are not known, 
a subroutine is used to compute approximate first harmonic flap- 
ping motions. If flight conditions but none of collective and 
cyclic pitch angles nor blade flapping motions are known, then 
collective and cyclic and blade flapping motions were- approxi- 
mated in preliminary calculations, using approximate uniform in- 
flow and such that thrust and drag trim resulted. 

If blade motions are not known, but collective and cyclic 
pitch angles and flight conditions are known, the blade motions 
used during the calculation of wake geometry are determined from 
the solution of the equation of motion for a rigid blade re- 
strained from flapping by a spring. The spring stiffness is zero 
for a rotor with a flap hinge; for a rigid rotor the spring stiff- 
ness should be adjusted to give a rigid-body natural frequency 
approximately equal to the first flapwise bending natural frequen- 
cy. Blade flapping motions are assumed to be forced by aerodynamic 
loads which are computed by assuming uniform wake-induced inflow 
and using specified values of the shaft angle and the collective 
and cyclic pitch angles . It should be noted that approximations 
to only the steady and first harmonic flapping motions are deter- 
mined by the above process and used in the wake geometry calcula- 
tions, and that quantities such as desired thrust (which is used 
in preliminary calculations) may not be the same as the computed 
thrust from this program. This is a result of blade motion effects 
due to blade flexibilities and nonuniform wake-induced velocities. 


The equation of motion of a spring-restrained rigid blade, 
may be written as 


d 2 g 


dt 


+ [k + ft 2 (I b + ) ] 6 = M 


( 1 ) 


where 1^ is the mass flapping moment of inertia about the hinge, 

t is time, 

m b is the blade mass, 

6 is the distance from the shaft axis to the hinge, 

x, is the distance from the hinge to the blade center of 

mass , 

k is the spring stiffness (zero for a springless hinge) , 

M is the aerodynamic moment exerted about the hinge, 

3 is the flapping angle, and 

ft is the rotor rotational speed. 
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The moment M is given by 

R 

M = a (r - 6) dr 
r o 

where a is the lift per unit span, and is given by 

l = ircp (fir + sin fit) (a 0 [fir + sin fit] - w 0 

- rB - V_B cos fit + V,a + ficB/2} 
r f s 

Here r is the distance from the axis of rotation, 
r 0 is the radius to blade root, 

R is the rotor radius, 
c is the blade chord, 
p is the air density, 

Vf is the forward flight speed, 

a q is the sum of collective and cyclic pitch angles and blade 
twist angle, 

a is the shaft tilt, positive aft, and 
s 

w 0 is the uniform wake-induced inflow, and is given by 
w o = fiR(.5[(C 2 + y ^ ) ] * 5 - y 2 } * 5 

where C T = (aircraft weight) / (irpfi 2 R 4 ) 
and y = V^/fiR. 

By letting 6 be assumed to have the form 

B = Bo + sin fit + B 2 cos fit , 

substituting into the equation of motion, and discarding all but 
steady and first harmonic terms, there result three equations for 
Bq/ $i, and B 2 which are solved to define the blade rigid-body 

steady and first harmonic flapping positions. 

The natural frequency of free rigid-body flapping motion is 
given by 

to = { [k + fi 2 (1^ + m^ Bx^)]/!^} / 
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and for a rigid rotor, the desired spring stiffness, k, may be 
calculated from 

k = - ft 2 (1^ + m^Sx^) (2) 


where uj is the first flap bending frequency of the blade which is 
being modeled. (For a cantilevered rotor 6 may reasonably be 
assumed to be the distance from the shaft axis to the first sig- 
nificant flapping flexibility of the blade.) 

Determination of wake geometry .- The right-handed Cartesian 
coordinate system used in the formulation for the wake geometry 
calculations has the x-axis in the downstream direction (positive 
toward the rear of the helicopter, <Jj= 0) and the z-axis parallel 
to the rotor shaft of one rotor and positive up, as indicated in 
Figure 3 . The center of one rotor is located at the origin of the 
coordinate system. A second rotor, if one is used in the model, 
may be located arbitrarily within this system. The flapping mo- 
tions for each rotor's blade are determined independently as dis- 
cussed previously. Given the flight conditions and appropriate 
rotor and control parameters, wake geometry calculations are 
essentially those required to compute vortex element strengths, 
vortex induced flow, and motion of the vortex elements. 


The Biot-Savart law (see ref. 13, for example) is the basic 
relationship used for calculating vortex induced flows, and gives 
the fluid velocity j at a point located by the vector s_ due to 

a vortex with circulation r as 


q (s ) 

— — p 


_1_ 

4 7T 


r (£) SjXd£ 


s 


3 

1 


(3) 


where ^ 



s . 


In order to compute total vortex induced flow at a point, the 
integral is taken over all vortex elements in the flow. For a 
straight element with end points located at A and B the induced 
velocity, q , at C is given by 


q 


w 


ra <cos s a - oos V 


(4) 


where r is the strength of the element and 6^, o^, and d are as 

shown in Figure 4 . When computing the vortex induced flow at a 
point due to an adjoining vortex element, equation (4) becomes 
indeterminate. Crimi , in ref. 10, gives the vortex induced flow 
due to adjoining finite core radii vortex elements whose end 
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points are used to define an appropriate radius of curvature, S, 
(called script m in ref. 10) of a circular arc by the expression 


q 


s 


i 


1 

8ttS 



8S 


i- 


1 


tan 



+ 


1 

4 




fas * 1 ] 

1 

| 

+ r . 

In 

a7 tan X 

+ 4 




1 J 

— 

J 


(5) 


where r, a, and <j> are the vortex strengths, radii, and angle as 
indicated in Figure 5. Where it is necessary to compute the in- 
duced velocity at the edge of a wake where there are not vortex 
elements on either side of the point being considered, such as 
that at the juncture of the shed and trailing vortex elements at 
the tip (or root) as indicated at point A of Figure 2 , the arc 
radius is computed from the two nearest elements ' end points , and 
the induced velocity includes only the effect of the nearest vortex 
element. For example, if in Figure 5 represented the edge of 

a wake with no element corresponding to P^_ to the left, then S 

would be the same as in equation (5) but the term involving 

would not be included. The induced velocity of an element end 
point just shed from the rotor blade includes a component due to 
the blade (bound) circulations, so that blade downwash effects are 
more nearly approximated in the region close to the blades. 

The interaction of vortex elements is governed by the above 
equations except where a vortex end point falls inside the finite 
core of another element. In such cases the theoretically predict- 
ed induced velocity may be unrealistically high, and not repre- 
sentative of the true physical situation. In order that this un- 
realistic situation does not occur, the maximum induced velocity 
that is allowed is arbitrarily restricted to be less than a speci- 
fied percentage of the rotor tip speed. 

The blade circulations are calculated as follows. The veloci- 
ties V and U, normal and tangential to the plane normal to the 
shaft axis, are given approximately by 

V = V f (o g - a ) + w (6) 

U = fir + Vf sin i p (7) 

where a is the shaft axis angle with respect to the free stream, 
s positive aft, 

a is the forward tilt of the rotor plane with respect to 
^ the shaft axis due to flapping, and 
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w is the induced upwash due to the wake (w=0 at start-up, 
i.e., there is no uniform downwash at start-up). 

The momentum due to the average wake-induced upwash was cal- 
culated for some cases, and was within ±15% of the uniform down- 
wash momentum value, depending somewhat on the radial distribution 
assumed from the outboard load point to the blade tip. 

Using small-angle approximations, the lift per unit span is 
given by 


£ = | p U 2 c c £a [a 0 + V/U] (8) 

where c^ a is the lift-curve slope, a constant, oiq is the sum of 

collective and cyclic pitch and blade twist angle, and c is blade 
chord. The "Kutta-Joukowski law" states that (assuming V/U<<1) 

£ = p u r b (9) 

where is the blade circulation. Therefore 

= .5 c c^ ^ [a o (f2r + V f sin ip) + V f (a g - at ) + w] (10) 

For a set of rotor blades which have stepwise radial and azimuthal 
circulation variations, the above equations may be thought of as 
applying to each radial and azimuthal location independently. The 
wake-induced velocity on the blade, w, is made up of velocities due 
to known circulations in the wake and to unknown circulations at the 
blade, and may be written in the form 


W(r i'\> = W N (r i'V + l l °£j (r i^k )r(r £^j } 

3 


(ID 


where w N (r <l>^) is the induced velocity due to all known wake circu- 
lations , 

r (r , ip . ) is the blade circulation at r , ip . , and 

36 J x, 3 

a .(r^,tf/,) is an influence coefficient which, when multiplied 

by the circulation r(r.,^.) gives the induced 

3 

velocity of that element at r^,^. 

The summations over the indices £ and j indicate a summation over 
all radial sections of all blades at their respective azimuthal 
positions. Then a set of equations for all r's may be obtained, 
and is of the form 
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2 c i,a % l °2.jik T Si: 
^ 3 


+ a 0 (^ r i + v f sin l i J k ) + v ft a s “ 


Here is equivalent to tpj,) , and occurs on both sides of 

the equation. This equation is solved with a simple iterative pro- 
cedure. The procedure is as follows. The terms on the right-hand 
side of equation (12) which do not involve the F's are used as an 
initial estimate for the r^^'s. Then the r^'s are evaluated from 

the equation, and as each is evaluated, it replaces the appropriate 
(on the right-hand side). Once a complete set of F^'s is 

determined, it is compared to the previous set, and convergence 
is assumed when the sum of the squares of the differences of the 
sets of r's divided by the sum of the squares of the F's is less 
than an input control constant, approximately 0.00005. 

After the blade circulation values are determined for a given 
azimuthal location of the blades, vortex-induced velocities are 
computed at all vortex element end points in the wake. The blades 
are then advanced through an azimuthal increment, , and each vor- 
tex end point which is not attached to the blades is allowed to 
translate for the time period Ai p/Q, with the resultant velocity 
due to the free stream and induced velocities . The entire compu- 
tational process for the new rotor blade positions is then re- 
peated. That is, new wake -induced and free stream velocities at 
the blades are computed and used to determine new blade circulations. 
Then new vortex end point velocities are computed, etc. Each blade 
advance results in an additional set of shed and trailed vortices 
being added to the wake. The number of revolutions of wake re- 
tained for actual computational purposes is restricted by an input 
to the computer program. 

Wake flow and wake-induced velocity influence coefficient 
calculations . - Calculation of wake flow is done by using the same 
basic programming which computes flow at a point due to an arbi- 
trarily located vortex element, except that the point is speci- 
fied by input, and the flow is averaged and nondimensionalized for 
purposes of comparison with experimental data. Calculation of 
input for the blade loads part of the program is done by manipu- 
lating and properly subscripting numbers equivalent to the 
o^j(r^,^^) and similar numbers used in the computation of w K . 

These numbers are used in the solution of an equation in the blade 
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loads program which is approximately equivalent to equation (12). 
Both the wake flow calculations and the evaluation of the input 
to the blade loads program are computed only after a specified . 
number of revolutions have been done. Thus, the input (including 
both a and r type quantities as discussed above) from the wake 
geometry program to the blade loads program is based on approxi- 
mate (specified) blade motions and an approximately repetitive 
wake. The repetitive nature of wake-induced effects was deter- 
mined by visual checking of data to occur after approximately 
ftR/Vir rotor revolutions. 


Blade Loads and Response Model and Formulation 

The right-handed coordinate system used in the calculation of 
(periodic) blade circulations and blade response is located such 
that the z-axis is fixed to the shaft, directed upward, the x-axis 
is downstream, blade azimuth angle, , is measured with respect to 
the x-axis, and the distance radially outward from the axis of rota- 
tion on a given blade is denoted by r , as shown in Figure 6. 

Determination of aerodynamic loading .- The aerodynamic loading 
at a given radial and azimuthal station is derived from the flow 
experienced by the blade section, as sketched in Figure 7. The 
geometric incidence of the section with respect to the rotor plane 
is the sum of the rigid body pitch angle 0 q(^)/ blade twist O^r), 

and the torsional deflection ^(r,^). The velocity component 
tangential to the rotor plane, U, is given by 

U = ftr + cos a g sin ip (13) 

and that normal to the rotor plane, V, is given by 

V = - V f sin + h(r,i|/) - w (14) 

where h is the plunging velocity due to the response of, the blade 
including blade flexibilities, and w is the wake-induced downwash. 

As is indicated in Figure 7 , the foil section is replaced by 
a vortex distribution of strength y (r,C,^) along the chord. This 
distribution is adjusted to make the flow at the section tangent 
to the chord, which relates y to the w and yields the basic rela- 
tionships governing the aerodynamic loading. The function y may be 
written as the sum of two contributions: y 0 , which is the circula- 

tion required for flow tangency with the section subject to a chord- 
wise uniform downwash, and yj, which is the circulation required for 
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flow tangency for the section subject to a downwash which is vary- 
ing chordwise in such a way as to give no average downwash. Thus, 
the chordwise average of the downwash, w, is defined as 


w(r,<f») 


1 

2b 


f b 

wCr,c,4>) d? - 


-b 


The portion of the downwash which is variable over the chord, wi , 
is given by 


w 2 (r , £ ,ip) 


[w(r, - vr(r,f) ] 



+ <f>) • 


Now the flow over an airfoil subject to a chordwise uniform down- 
wash is exactly equivalent to the flow due to angle of attack. The 
distribution y 0 , or rather its integral, r 0 , which gives the lift, 


is that quantity which is measured in a wind tunnel , under static 
conditions. That is, the total circulation, r 0 , may be defined in 


terms of yq as 


r o 


I" T. «. 


-b 


or in terms of the lift per unit span, £, for incidence a and Mach 
number, , as 


r 0 = £ (a,M N )/pU 


where U is the resultant free-stream speed. In terms of the section 
lift coefficient, c , 


r 0 = U b c^a,^) 


( 15 ) 


where c^CajMj^) is an experimentally determined lift coefficient. 

The circulation resulting from chordwise variations in downwash is 
not readily measured, but thin airfoil theory allows evaluation of 
r x , the circulation which results from the downwash wj , as 


r b 


= - 2 


w x (r,S ,i|>) [ (b + S)/(b - 5)] 


1/2 


as 


( 16 ) 


-b 
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where 


r l = 


Y i d£ • 


Then the blade circulation is given by 

rlr i 'h> ’ b i (U ii + V ii )1/2 <««'“«) 


- 2 


b. 


1 W] (r^e,^) [ (b i + C)/(b i - 


5)1 1/2 d ? 


(17) 


- b i 


where i and j range over all radial and azimuthal stations, respec- 
tively, and where 


U. . = fir. + V- cos a sin ip . , 

i r s j 

V ij = ~ v f sin “ s “ w(r ± ,i(;j) , and 

a ± j = 0 0 ('I' j > + (r ± , 4,^ ) + 0 1 (r ± ) + tan -1 (V^/U^) 


The average wake-induced downwash effects are given by 


= I I o (r i ,^ j > r (r ,* ) (18) 

m n J 

where the function G mn ^ r -L^j^ ' when multiplied by / repre- 

sents the cumulative contribution to w of elements in the wake 
having strength r(r rijj )• The integral over the chordwise-variable 

downwash is found to be 



w 1 (r i ,c,4< j ) [ (b ± + C)/(b i - 5 ) 1 1/2 d? 

Tib . 

= - sr [0 o^ j + 1 ) + 4> ( r i'^ j + 1 ) ~ 


$ (r i^j_i )] (19) 
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Equations (17) are a set of N R x N A nonlinear equations in the 
N r x N a strengths r(r ± ,i|>j). (N R x N A is the total number of cal- 
culation points in the system. ) Determination of the circulation 
is done by an iterative technique as follows. Blade circulations 
from the wake geometry program are used to compute the wake-induced 

downwash w. and are also, used to compute aerodynamic forces, from 
a 3 

which blade response in terms of normal modes is determined. (The 
response part of the program is discussed later.) This allows deter- 
mination of h, 4> , and other quantities which are a result of blade 
response. The velocities and V^j , and angle of attack, a^, 

are then computed, and for each combination of i and j, c (a. . ,M. .) 

^ J 

is determined from an analytical approximation to measured lift 
coefficient values. Determination of these quantities allows a 
new set of circulations to be computed from equations (17) . The 
process is repeated, i.e., the new values of r(r^,i|ij) are used to 

obtain new estimates for r i , w, and a, which in turn allows calcu- 


lation of another set of r (r^ , <j> ^ ) ' s . The iteration continues until 

the values of the circulations for two successive iterations agree 
to within a prescribed amount. 


Once blade circulations have been obtained, the lift and mo- 
ment per unit span are readily calculated. Resolving the forces 
into components normal and tangential to the rotor plane, F z and 

F , respectively, results in the following expressions for F and 

X " 

F (with F positive in the normal drag direction) : 

X X 


F z = pur + p dt 


b r + 2c 


- b 2 /u 2 +v 2 - 

m 


- w) A. 


(w - w) /b 2 -^ 2 d£ 


( 20 ) 


-b 


+ c.b V /u 2 +V 2 , and 

d 


F = - pVr + p c, b U Xj 2 +v 2 
x d 


( 21 ) 


where c m and c d are the section moment coefficient (about midchord) 

and the section drag coefficient, respectively. The moment about 
midchord, M 0 , is given by 


M 0 


PU 


f b 

Y(?,t) 5 d? 


-b 


P _d_ 
2 dt 


(b 2 


? 2 ) Y (5 ,t) d£. (22) 
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The above quantities are functions of radial and azimuthal position, 
and are understood to have appropriate subscripts. The aerodynamic 
coefficients are computed from an empirical analytical representa- 
tion of experimental data for a NACA 0012 airfoil., and are non- 
linear functions of Mach number and angle of attack. Representa- 
tion of other airfoils may be implemented by developing subroutines 
which, given angle of attack and local velocity and Mach number, 
provide the necessary aerodynamic coefficients (lift, drag, and 
moment) . 


Determination of blade response . - The aerodynamic forces and 
moments , as computed, are not in a convenient form for calculation 
of blade response. A coordinate transformation and integration of 
the aerodynamic loads are performed so that generalized forces act- 
ing on previously determined flexible blade normal mode shapes may 
be computed. The moment is transferred to quarter-chord from mid- 
chord and the nondimensional distributed airloads are made dimen- 
sional. . The distribution of aerodynamic forces and moments is 
approximated by straight lines passing through zero at the blade 
root and tip, and through the computed values at the aerodynamic 
loads points between the root and tip. Lumped loads to be applied 
the lumped parameter blade model mass points are then computed 
by integrating the distributed loads over appropriate radial lengths. 

4- Tn 

For example, the drag force applied to the i — mass is given by 


f 


x 


■R2 

F^(r) dr 
r i 


where R^ and R^ +1 are midway between masses nu and itn , and between 
nn and m ±+i , respectively, and where F^(r) is the straight line 
approximation to F x . Loads in the local blade coordinate system 
result from the coordinate transformations 


F =-f cos $ - f sin $ 

V z X 

F = f sin $ + f cos $ 
w z x 

M, = m n + z„ F 

<p u A V 

where z^ = distance of elastic axis forward of quarter chord. . The 

total forces and moments at the radial positions corresponding to 
each blade mass point and azimuth are computed. (These loads are 
applied at the blade elastic axis.) For example, the force at 
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the i th. mass station and jth azimuth., in the direction of blade 
flapping motions, denoted by Q , is given by 

V ij 

Q v _ = iru(fi 2 sin (h^ + e^ cos $^) 


e c U i + h i 


cos 


($ . . 

ID 


e 0 ) /cos $ ± .J 

j 


- 2Q,z . 
1 


*iD 


sin 




+ 



(23) 


Here the first term is the steady acceleration force on the i th 
mass due to rotation, the second term accounts for control angle 
acceleration, the third term is a gyroscopic coupling term, and the 
fourth term is the aerodynamic load on the ith mass. Equivalent 
expressions are used for Q^, Q^, and , which correspond to the 

elastic deflections V , w, and <|> , the in-plane angular, in-plane 
linear, and torsional angular elastic deflections, respectively. 

The Q is given by 


Q ?ij - »i ° 2 E i r i 


- 2 ft (I 0> sin + itk sin B 0 /cos $ i j)0 c 


- 2ft [I o 5> i; . - m ± (v ± j sin cos )], 


(24) 


Q w is given by 


w 


ID 


nu ft 2 cos (h^ + cos 


- [m. h. sin ($ .. - 6 0 )/cos 4> . . ] 6 

IX lj j XJ L. 


- 2 m.fie. cos 4> . . V . . + F 
li i] id w 


ID 


and 


(25) 
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Q. is given by 

Q = - fl 2 sin <&..(I 0 cos + m. e. h. ) 

+ ij 13 i 13 iii 

- [i 0 _ + m i e ± h ± COS - e 0 )/cos $ i j]e c 


+ 2 fi I 0 sin $ . . v . . + M 
13 ID 


♦i 


( 26 ) 


In equations (23) through (26) 


m^ = mass at i th radial station. 


n = I + m. e. 

i x i 


= torsional inertia about c.g. at ith radial station. 


e . = chordwise separation of the elastic axis and center of 

1 mass. 


r. = radial distance to the ith mass, 

1 — 

= angle between the chord and the rotor plane, 

h. = horizontal separation of the elastic axis and the 
1 pitch axis, 

0 0 = blade angle at blade root, 

j 


0 = pitch input control angle, and 

w, v, Y, and <j> are values from the previous iteration. 


The generalized force associated with each mode, F R , is then com- 
puted according to 


F K (t) - 


JTo a + Q A + Q a + Q, A, ~] + l o v 'r, v ^ 

iL~ V i V iK W i W iK *iK +i *i K J j KJ K1 


(27) 


where the A are mode shape deflections in the v direction, at 
v iK 
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the ith mass station, for the Kth mode, etc., and the a 


KJ ' K' 


and 


5 are the artifical damping coefficient, natural frequency and 
generalized coordinate rate for the previous iteration, respec- 
tively, for the K th mode. (The mode shape quantities, A , etc., 

v iK 

and corresponding natural frequency are provided by an auxiliary 

program. The corresponding model and formulation is briefly des- 
cribed later.) The generalized coordinate in the Kth mode, 

and its time derivatives are then computed, based on the equation of 
motion 


C K + 2 °K W K + W K ? K = F K^ 


(28) 


and the requirement that £ K be periodic. Response variables of 

interest, such as the flapping velocity due to elastic motions, 
V^(t), are then computed from expressions of the form 


v ± (t) = l 


K 


iK 


t K (t) 


(29) 


Blade dynamic response quantities are then used to define, by appro- 
priate transformations and interpolations, quantities to be used 
in blade load calculations such as h(r^, \p . ) , as given by 


h = {v cos $ - w sin $ + z^ <}> cos $ 


+ [z a - h][e - fi(£cos 0 + <p sin <J>) ] } 
A C 


+ y cos a T cos at (£ cos $ + \p sin <f) 


(30) 


The process is then repeated, with new blade circulations 
being computed with the inclusion of the effects of the new blade 
motions. The iteration process is terminated when (1) a specified 
maximum number of iterations between blade circulations and blade 
motions have been performed, (2) the values of blade circulations 
for successive iterations agree to within a certain constant, (3) 
the iteration procedure begins to diverge, or (4) the time or output 
line or page limit as controlled by the computer program job card 
are exceeded. Upon convergence, divergence, or reaching the limit 
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on the number of iterations of the iteration procedure between . 
blade circulations and blade response, shear and moment quantities 
are computed from mode shape quantities. For example, the torque, 
T, is given by 


T (r ± , ^ j ) 





(31) 


Both aerodynamic forces and shear and moment quantities are harmon- 
ically analyzed, and are printed in terms of harmonic coefficients 
as well as functions of radial and azimuthal position. 

Determination of initial estimates .- Due to the iterative solu- 
tion for a compatible set of blade loads and blade response, good 
estimates of load and response quantities are helpful to prevent 
divergence, and if properly chosen enhance convergence. 

Estimates of blade flapping motion for use in blade loads 
calculations may be computed by a method equivalent to that used 
in determining blade flapping motions for the wake geometry and 
wake-induced velocity influence coefficient calculation. When 
the shaft axis is chosen to be perpendicular to the tip path plane, 
blade flapping motions are assumed to be zero. Estimates for 
blade (bound) circulations for the blade loads program are input 
from the wake program, and are those values associated with the 
last steps of the wake program calculations. 

Blade natural frequency and mode shape determination model 
and formulation . - The program whichTcalculates blade natural fre- 
quencies and mode shapes is an independent program which was de- 
veloped by RASA and could be replaced by any approximately equiva- 
lent program. The necessary output for any such program in order 
that it might be used as input for the blade loads program are the 
natural frequencies together with the corresponding mode shape 
quantities, i.e., flapping and edgewise displacements, angles, 
shear forces, and moments, and torsional deflection and torque. 

These mode shape quantities are to be defined at the location of 
the point masses of the lumped parameter model, with mode shape 
magnitudes adjusted to give unity generalized mass. The lumped 
parameter model lengths, mass and inertias, mass eccentricities, 
offsets of elastic axis from pitch axis and midchord, twist distri- 
bution (but not model bending or torsional stiffnesses) are used 
in the blade loads-response program. 


The blade frequency program developed by RASA is described 
in references 14 and 15, and the model used, development of input 
data, and operation of the blade frequency program are discussed 
in more detail in those reports. Briefly, the model used for the 
real blade is a lumped parameter approximation consisting of 
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uniform massless elastic beam sections under tension due to centri- 
fugal loads, with point masses and inertias located at the ends of 
the massless lengths. A modified transfer matrix approach is used 
in determining the natural frequencies and mode shapes. (See, for 
example, ref. 16) . The natural frequencies and corresponding mode 
shape quantities are used in the calculation of the response of the 
flexible blades to aerodynamic and inertia loads in the blade loads 
and response part of the program. 


Computer Program Description and Use 

The computer programs used for the implementation of the model 
constitute more than one separate program, and may be classified 
as those used for preliminary calculations, that used for wake 
geometry and wake flow or wake influence coefficient calculations, 
and that used for computation of blade loads and blade response. 

The preliminary calculations, some of which were done by computer 
programs, included determination of blade control inputs and deter- 
mination of blade dynamic properties as described by a limited 
number of fully coupled free vibration mode shapes and correspond- 
ing natural frequencies. The flow diagram shown in Figure 1 indi- 
cates the general procedure used in the calculation of results. 
After available flight conditions, control settings, blade proper- 
ties, etc. are obtained, this data is used to develop input to the 
wake program and the blade loads and response program. Equations 
which were programmed for calculation are given in the Appendix. 

The necessary data for input, including program control vari- 
ables, is given in the Program Documentation, which is available as 
RASA Report 69-9. Experimentally measured or theoretically or 
empirically calculated data may be used for such data as blade 
cyclic and collective pitch settings, blade flapping motions, blade 
natural frequencies, etc. The blade control settings and approxi- 
mate flapping motions for cases where such data is available may 
be taken from such data. Where experimental data isn't available, 
performance-type programs, charts, or tables may be used to pro- 
vide this data, as appropriate. Blade mode shapes and frequen- 
cies may be computed by any method desired, but must give mode 
shapes with a generalized mass of unity for use in the blade re- 
sponse calculations. The blade vibration characteristics used in 
the calculations discussed herein were obtained with an independent 
program which computes fully coupled (flapwise, edgewide, torsion) 
free vibration mode shapes and natural frequencies. Airfoil 
characteristics used in the calculations were for a NACA 0012 air- 
foil, and are determined in a subroutine by algebraic approxima- 
tions to experimental data. If desired, this subroutine could be 
replaced by similar subroutines or table look-up procedures for 
other airfoil sections, for models which include rate effects, 
etc. The calculation or tabulation of airfoil data in a form 
suitable for use in the wake and loads programs may be thought of 
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as a preliminary calculation, since at the present time it is 
operational for a NACA 0012 airfoil section only , and any differ- 
ent models would have to be developed in a similar manner before 
being used. 

After the necessary input quantities have been determined, the 
wake geometry program may be run to obtain wake geometry data and 
wake flow data or wake influence coefficients for use in blade loads. 
Ordinarily, if wake flow data is desired, wake influence coeffici- 
ents are not, and vice versa, but both may be obtained if desired. 

The blade loads and response program does not require wake 
influence coefficients for operation; omission results in blade 
loads and response which do not include free wake effects; rather 
an assumed uniform wake-induced inflow is used. It is necessary 
to input at least one natural frequency and mode shape from the 
blade free vibration analysis, as well as flight conditions, heli- 
copter blade control settings, and program control variables. The 
blade loads and response program may be conveniently thought of as 
two large subroutines , one which computes blade loads , and one 
which computes the response to those loads. Nonlinear effects are 
included in the model, so starting estimates for some quantities 
are made either by input or by the programs , and iterations between 
load and response calculations are required in order to compute a 
compatible set of blade loads and response variables . The maximum 
number of iterations allowed is limited by input. Convergence 
within that number of iterations is determined by the relative 
change in bound circulation from one iteration to the next, as com- 
pared to an input number. Output of load and response data is pro- 
vided when any of the following situations occurs: (1) the conver- 

gence criteria between loads and response is satisfied, (2) diver- 
gence of loads from one iteration to the next occurs, or (3) the 
maximum number of iterations is reached. This output data in the 
first situation is the desired final output, while in the other 
cases the output serves as a source of information which may be 
used to determine reasons for failure to converge. 

Listings of the program used to compute the wake geometry and 
wake flow or wake induced velocity influence coefficients, and the 
program used to compute blade loads and response are given in 
Volume II. Listings of the programs used to perform preliminary 
calculations, including the blade free vibration program are not 
presented, inasmuch as they were not part of the contract effort 
and are proprietary programs . It should be noted that the program 
input and output requires the use of English unit quantities, 
rather than metric unit quantities. 
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RESULTS AND DISCUSSION 


This section presents and discusses the results of prelimin- 
ary calculations, the comparison of predicted and measured wake 
flows and the effects of various numerical parameters of the pro- 
gram on the representation achieved by the model.. The predicted 
blade load and response quantities for a model rotor are then pre- 
sented and discussed. The effect on blade loads of geometric and 
flight conditions including advance ratio, vertical blade separa- 
tion, azimuthal blade separation, and blade radius ratios are dis- 
cussed. All calculations of blade load and response quantities 
include the effects of a free wake. The free wake and correspond- 
ing approximate blade loads were computed with the same basic pro- 
gram which was used to compute wake flows. Blade loads and re- 
sponse for flexible blades as computed included the effect of a 
free wake for specified blade motions, and the effects of flexible 
blades as contained in the mode shapes and natural frequencies of 
a fully coupled (flapwise, torsion, edgewise degrees of freedom) 
lumped parameter model of a rotating beam. Various wake parameters 
chosen during and prior to the wake flow runs (which were made for 
comparison of predicted and measured flows) , such as number and 
spacing of trailing vortices and . azimuthal step size in the full 
mesh and modified wake, were used for subsequent wake geometry 
and blade loads calculations . The problems involved in the use 
of a full mesh model for the entire wake have been mentioned 
previously, and will be discussed in more detail hereafter. The 
blade free vibration calculations include the effects of rotor 
speed (flapwise and chordwise centrifugal stiffening) , nonuniform 
mass, inertia, twist, and stiffness distributions, and spanwise 
variations of the center of gravity and the elastic axes with re- 
spect to the pitch axis. It does not include the effects of steady 
blade lag or coning (i.e., the free vibration equilibrium position 
has zero lag and coning angles):, nor gyroscopic coupling terms 
(which would necessitate the use of complex numbers in the calcula- 
tions, and would result in complex mode shapes) . Gyroscopic coup- 
ling terms were accounted for in the calculation by including them 
in the generalized forcing function for each mode in the blade 
loads - response portion of the computer program. 


Results of Preliminary Calculations 

Preliminary calculations for the wake geometry and blade loads 
calculations included determination of a lumped parameter model of 
the various rotor blades and the corresponding fully coupled natural 
frequencies and normal modes. These calculations were performed by 
an independent program, and results were used as input in the blade 
loads-blade response program. The various blade models are sketched 
in Figure 8, blade properties are given in Table 1, and their 
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corresponding natural rrequencies and mode shapes are given in 
Table 2 through 5. The changes in blade models are based on 
changes in tip radius or airfoil cutout. The first spanwise 
section torsional stiffness value, includes the effects of a con- 
trol stiffness adjusted. so as to. give a control torsion frequency 
ratio, oj/fi, of 3,50, for the standard blade. The second edgewise 
bending mode is in the frequency range of the other modes only for 
the case of the blade which has an airfoil cutout to 0..4R. However, 
for completeness this mode was calculated and used. in the early 
blade response computations for all blades. It was evident in the 
first calculations for case 9, that the relatively low chordwise 
and torsional bending stiffnesses of . blade model 4 inboard of 
r=0.4R (given in Table 1) resulted. in excessive torsional and 
lead-lag deflections. The intent of the calculations with this 
configuration was to determine the aerodynamic effects of such a 
rotor configuration. Thus, it was. decided that the blade stiff- 
nesses , as contained in the blade normal mode shapes and frequen- 
cies, should be those of the blade with airfoil cutout to r=0.1, 
so that aerodynamic effects, and not elastic or dynamic effects 
would dominate. The "very flexible" blade properties as given for 
case 9 in Tables 1 and 5 are given for completeness, but were not 
used in the final blade loads, rufts of this configuration; rather, 
mode shapes and frequencies as determined for blade model 1 were 
used, together with airfoil chord size, cutout to r=0.4, etc., 
from blade model 4 . 

A performance-type calculation was made to determine collective 
and cyclic pitch and first flap bending amplitude for all of the 
flight conditions except case 1. Experimentally measured data was 
used to obtain input parameters for case 1. The model used in the 
performance calculation is essentially an isolated rotor for which 
values of collective and angle-of-attack of the plane of no feath- 
ering were adjusted to give approximate drag and thrust trim, 
with fuselage drag included. Cyclic pitch angles in the shaft 
oriented axis system were determined from the angle of attack of 
the plane of no feathering and an assumed shaft angle. The shaft 
was assumed to be perpendicular to the tip-path plane in subse- 
quent calculations. Quantities used as input for the free wake 
and blade loads calculations are given for various flight condi- 
tions in Table 6. These include blade model identification (as 
given in Table 1 and Figure 8) advance ratio, blade chord, vertical 
separation (Az) , azimuthal separation ( A ) , '.shaft angle, and collec- 
tive and cyclic pitch control angles. Two rotors were used in all 
calculations, and these were numbered 1 and 2. Rotor number 1 has 
the standard blade tip radius in all cases, and Az and A are posi- 
tive for cases where rotor number 2.v. is positioned above and ahead 
of rotor number 1, respectively. Other quantities are as defined 
in the previously discussed figures and equations. For case number 
11, the collective and cyclic pitch angles are given for rotors 
number 1 and 2, in that order. All other cases were run with the 
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same collective and cyclic values for both rotors. (It may be noted 
that cases 1 and 2 could have been run as a 4-bladed single rotor 
instead of as 2-bladed double rotors.) All wake geometry and 
blade loads calculations were run for NA=16 in which the blade 
azimuth interval is 22.5°. It was found during the first blade 
loads run of case 7 that the rotor was badly out of trim. The 
performance calculations for case 7 were greatly affected by 
stall, which resulted in very high computed values of collective 
and lateral cyclic pitch. In contrast, the blade loads calcula- 
tions indicated considerably less stall, revealing that lower 
values of collective and lateral cyclic pitch were necessary. 

These lower values were then used in the final run for case 7 , 
and are those given in Table 6 . Downwash momentum values were 
approximately the same as theoretical values , depending somewhat 
upon assumed distribution from the outboard load point to the 
blade tip. 


Full Mesh Wake Model Results 

Wake flow with fu'll mesh model .- As mentioned previously, 
predicted and measured wake flow comparisons revealed such serious 
undesirable characteristics of the full mesh model that a change 
in the wake model was made. The nature of the discrepancy between 
the predicted and measured wake flows can be seen in the plot of 
the time-averaged wake induced flow shown in Figure 9. This shows 
unacceptably large variations in the predicted flow in the inner 
portion of the rotor disc along a port to starboard line at x=0 
(i.e. , at the rotor hub). No such variations. were observed at 
0.5R forward of the rotor hub, but similar variations were pre- 
dicted along a line at 0.5R aft of the rotor hub. After detailed 
investigations of both the computer program and output, it was 
concluded that the large fluctuations in predicted average wake 
flow were due to shed vortex elements in the wake model. The shed 
vortex elements which contributed to the fluctuations were approxi- 
mately aligned in the free stream direction and, because of azi- 
mutal variations in bound circulation, had sufficient strength to 
cause the significant velocities. For example, shed elements 
contributing to large fluctuations at position x=0.0R and x=0.5R 
were noted, in general, to be those shed from the blade as it 
moved in the region t|i=135° - 225° which in turn were carxied back 
by the free-stream velocity to these positions. The alignment of 
the vortex elements in the free stream direction allows the ele- 
ment to induce velocities in the same direction for the entire 
time over which the average is taken. Two methods for alleviation 
of this unrealistic characteristic were considered.. The full mesh 
wake model would be a satisfactory model provided that a suffi- 
ciently small azimuthal step size were used to yield shed vortices 
of low strengths and corresponding low induced velocities. The 
number of shed vortices needed to accomplish this was estimated 
to be greater than fifty, and the running time (or costs) of such 
a model would be prohibitive. Previous investigations have indi- 
cated that the shed vorticity immediately behind the blade is 
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important in computation of blade loads on that blade which sheds 
the vorticity, but is not as important in computation of loads on 
other blades. Dropping the shed vortices and reducing the number 
of trailing vortices was the second alternative which was investi- 
gated. The tip (trailing) vortex, of course, is known to be im- 
portant for computation of all blade's loads wherever blade- 
vortex interactions occur. Previous investigations indicate that 
inboard wake representation can also be important. 


Modified Wake Model and Results 


The model, as shown in Figure 2, which was successfully used 
to overcome the problem of large fluctuations in average velocity, 
was one in which the full mesh wake was retained for a few azi- 
muthal steps behind each blade, and then replaced by a modified 
wake, containing only trailing vortices. The number and location 
of trailing vortices in both the full mesh wake and in the modified 
wake are arbitrary, within dimension restrictions, and are deter- 
mined by input. Other changes in the wake model include program 
modifications to account for vortex stretching effects, vortex 
element representations of vortex sheets, and vortex element in- 
duced velocity limits for wake flow and blade loads effects. 

Vortex stretching . -The vortex stretching effect is that asso- 
ciated with a vortex which has been deposited in the wake with a 
length approximately rAi|>, and which changes length due to free 
stream and wake-induced velocity effects. For high advance ratios, 
the inboard vortex elements on the advancing side tend to get 
stretched significantly as they are deposited in the wake, and 
such stretching results in stronger induced velocities from these 
elements than experiments indicate exist. The program now accounts 
for vortex stretching by multiplying induced velocities from in- 
dividual elements by the ratio (&o/£) where q ± s the length rAt|> 

(the approximate length of a vortex trailed from radial position r, 
for azimuthal step sizes of Ai|0 and l is the current element length. 
This modification showed some beneficial effect at |i = 0.235, and 
is expected to provide increasingly more realistic wake models with 
increasing advance ratio. The model now approximately preserves 
the total vorticity associated with each element, where before the 
vorticity per unit length was conserved. 

Vortex core-smearing . -Use of vortex elements for the inboard 
wake gave correct trends in wake-induced flow, but resulted in 
peaks where such vortices were transported by a wake flow calcula- 
tion point in an approximately fore-and-aft orientation. This 
effect is observed experimentally for the tip vortex, but not for 
the inboard wake. A comparison of the induced velocity field from 
a set of vortex sheets with that computed for a set of finite core 
elements having the same total strength indicated that a realistic 
finite core vortex element representation of a vortex sheet was 
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obtained by using vortex core radii with the dimensions of one- 
half the sheet span which they represent. The wake-induced veloc- 
ity field for a set of four sheets of uniform strength and four 
finite core vortices with the same total vorticity is shown in 
Figure 10, and demonstrates the improved representation of a sheet 
by finite core vortices which may be expected with this model. 

The induced velocities as shown in Figure 10 are for a distance 
from the plane of the sheet of one-fourtieth of the sheet width. 

For distances larger than this , there is improved agreement between 
small core element and sheet-induced velocities, and at a distance 
of one-fourth the sheet width, the two velocity plots are coinci- 
dent. At distances less than one-fourtieth the sheet width the 
disagreement between smeared core and sheet-induced velocities be- 
comes more pronounced near the edges of the sheet. The agreement 
between sheet and smeared-core velocities in the interior remains 
good for essentially all distances from the sheet. In general, it 
may be seen that the usual finite core vortex elements give two 
large peaks of opposite sign near the center of the span which it 
represents, while the smeared core and sheet both give peak veloci- 
ties near the sheet edges. Differences between smeared core and 
sheet velocities in the interior region tend to be reduced by the 
cancellation of the large velocities due to neighboring sections 
having these velocities with opposing signs at the section inter- 
face. 


Inboard modified wake vortex core radii are determined by in- 
put so that this or other models may be used as desired, provided 
that the radii do not have azimuthal variation. Wake flow runs 
showed some improvement with this modification, and it is expected 
that the improvement would be more significant for flight conditions 
such that the inboard wake passes close to the blade load points, 
wake flow points, or other wake element end points, where the in- 
duced velocity from the "smeared" core vortex element model will 
be a better representation of the inboard sheet than the previously 
used "small" core vortex element model. (The tip vortex core has a 
"small" core, and is determined as described in the Appendix.) 

Induced velocity limit . -The use of a velocity limit for the in- 
duced velocity from an individual vortex element is a feature which 
was found to be necessary in order to allow computation of reason- 
able average velocities at wake flow points, and which is also ex- 
pected to give reasonable blade loads. (It is suggested that this 
limiting should be done differently for future runs where blade- 
vortex interactions with the blade and vortex parallel or more 
nearly parallel than perpendicular are expected to occur.) Figures 
11 and 12 show velocity time histories with and without velocity- 
limiting which, due to the finite azimuthal step size, result in 
different averaged velocities. The non-limited-velocity curve, as 
sketched, is more realistic, but the averaged velocity, v is biased 
by the location of the vortex element with respect to the flow com- 
putation point near ip = 150°, i.e., v is incorrect due to the time 
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spacing used to obtain v(t) values to compute v. The limited v(t) 
curve, while not physically realistic, gives a more acceptable 
value of average velocity. The same general behavior is expected 
to be shown for blade loads. That is, values of wake-induced 
velocities are determined only a finite number of times, at equal 
time increments. Use of a small number of times (or azimuthal 
steps) is essential for economical program operation, but may re- 
sult in an incorrect., average velocity, an incorrect loading for 
blade loads, and incorrect wake distortion. The type of inac- 
curacies which occur with computation of v are seen from Figure 
11. The blade loads inaccuracy, with the use of normal modes 
as generalized coordinates, is one whereby the use of a small 
number of azimuthal steps may tend to cause lower frequency modes 
to be excited more than would be excited by the real vortex-blade 
interaction. Since this study did not include the use of high 
frequency modes, the same type of vortex-induced velocity limit- 
ing as was used to calculate v is appropriate to keep lower modes 
from being incorrectly excited in the blade loads program. The 
wake flow and blade loads physical situations where this limiting 
is necessary are where a vortex passes by a wake flow or blade 
load calculation point such that the vortex is approximately per- 
pendicular to the free stream or parallel to the blade, respec- 
tively, and where the wake flow or blade load calculation points 
are near the vortex core at one calculation time. This situation 
was present for wake flow calculations on the advancing side just 
below the rotor disc. Use of velocity limiting gave much more 
realistic v versus space plots for those situations where vortices 
passed close by wake flow points, and did so for a reasonably 
small number of azimuthal points. 

The magnitude of the velocity limit is arbitrary as far as pro- 
gram operation is concerned, and is input as a function of the 
radial position of the vortex element. In the wake flow runs the 
velocity limit used was five times the momentum value of uniform 
wake- induced velocity, Vq. This velocity is a function of thrust 
coefficient, advance ratio, and rotor control angle, and is defined 
(see ref. 17) as 


V 0 /Sm = 1/2 C T /(1 _ 3/2y 2 ) (p 2 + X 2 ) ^ 


It is felt that this is a reasonable limit, based on the wake flow 
calculation results. It is suggested that a multiple of Vq be used 
for both wake flow and blade loads types of runs. (A limit of 10 
v 0 was used in the blade loads runs.) 

The computer program used to obtain the results reported here 
allowed considerable flexibility in choice of various model charac- 
teristics. It should be noted that the same model was used to ob- 
tain wake flow results for three different advance ratios, and 
that only parameters which depended on advance ratio were changed 
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for different advance ratios. The following assumptions were used 
in the wake flow and blade loads calculations , and were determined 
to be reasonable on the basis of both theoretical and experimental 
data. 


Model parameter determination .- Wake flows and wake-induced 
effects for blade loads were assumed to be periodic by the time 
that wake elements deposited at the forward edge of the rotor disc 
reached the aft edge of the rotor disc. In terms of the number of 
revolutions of the rotors, N, this occurs for N = 1/yir. The 
velocity-limit of 5 Vo was observed to give satisfactory results 
for wake flow calculations . Input values were not varied with 
advance ratio, except for variables which depend on advance ratio, 
such as thrust coefficient, cyclic and collective pitch, etc. 
Inboard vortex element core sizes were "smeared" in the modified 
wake, vortex stretching effects were included, the same number 
and location of vortex elements in the full mesh and modified 
wake, and the same number of azimuthal steps were done for all 
the advance ratios. Various combinations of these parameters 
were used in preliminary wake flow runs, and the effects of these 
variations are discussed later. Velocity contributions from dis- 
tant vortices were neglected in the wake geometry runs used to 
generate input for the blade loads runs. Preliminary calculations 
showed that such approximations were both reasonable physically 
and advantageous economically. 


Various methods for determining the location and number of 
trailing vortices in the full mesh and modified wakes were used. 

The final method for determining full mesh trailing vortex loca- 
tion was one which assumes a loading and velocity distribution 
linear in radial coordinate, divides the blade into NR sections 
with equal total load, and locates the trailing vortices within 
each section so as to conserve impulse of vorticity. The modified 
wake trailing vortices were found to be best located by placing 
the tip vortex at the full mesh tip vortex, replacing the full 
mesh root vortex by two modified wake vortices located so as to 
conserve impulse, and placing one (or more) additional vortices 
so as to conserve impulse of the remaining sections of the full 
mesh wake. The number of full mesh and modified wake trailing 
vortices is arbitrary, as specified by input; but present dimen- 
sions allow 5 or fewer full mesh vortices and 4 or less modified 
wake vortices. These numbers of vortices were used for the wake 
flow runs, and appear to give satisfactory results. Fewer trail- 
ing vortices in the full mesh would make radial blade load varia- 
tions difficult to predict, while fewer trailing vortices in the 
modified wake could result in unrealistic wake flow in the interior 
region immediately below the rotor disc. 
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Wake Flow Results 


Wake- induced flow in the region immediately below the rotor 
disc was computed and compared with the experimentally measured 
average flow as reported in reference 17. The plots are for two 
vertical positions for each of the three survey lines shown in 
Figure 13, as indicated in the figures, and present the nondimen- 
sionalized average velocity as a function of y, y being the star- 
board to port coordinate, +y being on the advancing side of the 
rotor, and all dimensions being nondimensionalized by rotor radius. 
Computations for advance ratios of 0.097, 0.141 and 0.235 were 
performed and the resulting plots of average wake-induced velocity 
are shown in Figures 14 through 22. Calculated velocities com- 
puted for a single tip vortex as reported in reference 10 are 
shown for the 0.141 advance ratio case. The overall agreement be- 
tween computed and measured average wake-induced flow for the 
present model is good. Both the region near the disc edges and 
the inboard region have now been quite accurately modeled, with 
the inclusion of a more realistic inboard wake representation. 

The improvement achieved by including the inboard wake is particu- 
larly noticeable when compared with the tip-vortex-only model re- 
sults at x = 0.5 and z = - 0.07, as shown in Figure 19. 

There is generally better agreement between computed and mea- 
sured flow along the forward survey lines than along the mid-rotor 
and aft survey lines. The presence of the rotor shaft and of the 
other structures associated with the experimental tests was not 
modeled in the calculated flow. The presence of the rotor mast 
precluded measurements near y = 0 at the mid-rotor survey line, 
and could influence the apparent disagreement between calculated 
and measured flows in this region, which disagreement appears to 
be worst for y = 0.235. The largest discrepancies between computed 
and measured flows are noted along the aft survey line, particu- 
larly for y = 0.097 and y = 0.235. In both cases the flow magni- 
tude is underpredicted near y = - 1.0, and for y = 0.097 there is 
a predicted "bump" between y = - 0.5 and y = - 0.7 which has the 
same trend as the measured flow, but which is much higher. The 
wake is highly distorted for this case, and it is possible that 
some of the wake distortion effects are not adequately represented, 
leading to the prediction of this "bump" . Use of more azimuthal 
steps with corresponding changes in averaged velocity values may 
more adequately represent the flow in this case. For the y = 

0.235 case, the disagreement near y = - 1 could be due to incor- 
rect wake flow point location. The measured data was reported to 
be taken parallel the tip path plane, but lateral rigid body flap 
was measured. The predicted flows were computed by assuming no 
lateral flapping and placing the wake flow points parallel the 
resulting tip path plane. If, in fact, the measured data was 
taken in a plane not parallel the- tip- path plane, the small rela- 
tive lateral flapping could have resulted in flow asymmetries that 
are observed in the measured data at y = 0.235, which are much 
more severe than for the other advance ratios . 
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Velocity limiting effects .- The effect of velocity-limiting 
for wake flow comparisons may be seen by comparing the predicted 
flow shown in Figure 19 (for x = 0.5, z = - 0.07) with the pre- 
dicted flow shown in Figure 23 . On the advancing side , the 
rather severe variations in v/v 0 with y have been traced to the 
tip vortex-induced velocity being close to some of the wake flow 
calculation points for only one azimuthal step,. The resulting 
v/v 0 quantities, as discussed in the previous section and shown 
in Figures 11 and 12, result in the variations shown in Figure 
23, and were removed by limiting the induced velocity contribu- 
tion at any computational point to have magnitude less than or 
equal to 5 v 0 (v 0 is the momentum value of uniform inflow) , in all 
cases. This limit, though somewhat arbitrary, was used for all 
flight conditions, and appears to be a realistic limit. 

The wake-induced flows shown in Figure 23 do not contain 
other model characteristics which were used for the final wake 
flow calculations, such as vortex stretching and vortex smearing. 

The effects of the variation of various parameters on wake flow 
calculations were noted during preliminary wake flow runs. It was 
observed that the number of azimuthal steps influences the loca- 
tion of v/vg variations due to vortex location when velocity limit- 
ing is not done, and that a significant increase in this number 
(from 12 to on the order of 100) is necessary to obtain realistic 
v values without velocity limiting when the wake flow points are 
such that a tip vortex passes close by them at an oblique angle. 
Running time is increased significantly with an increased number 
of azimuthal steps. The number of radial trailing vortex elements 
required to compute reasonable radial aerodynamic load distributions 
is felt to be five or more, five being the upper limit of the pres- 
ent program due to core storage restrictions. The number of trail- 
ing vortices in the modified wake was varied from one to four, 
and the placement of these vortices was also varied. One trailing 
vortex gives the best results when placed near the tip, but does 
not consistently predict the flow in the interior of the disc 
immediately below the rotor, as can be seen from Figures 17, 18, 
and 19. Use of two trailing vortices generally resulted in poor 
prediction of the interior flow also, with either very little 
difference from the tip vortex only model, for small circulation 
values of the inboard element, or large spacial variations in v 
for large circulation values. Three trailing vortices gave similar 
results. It should be noted, however, that runs with only two or 
three trailing vortices were made without vortex smearing, and 
vortex smearing could improve the results considerably for such 
cases. The best representation was obtained by using four trail- 
ing vortices in the modified wake model previously described. 

Other models, in which the full mesh root vortex was not divided 
into two modified wake vortices, gave large spacial variations in 
v near the location where the root vortex is in an approximately 
fore and aft orientation. Again, no results were obtained 
with only vortex smearing but without division of the full mesh root 
vortex . 
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Calculated Wake Geometries 


Free wake effects .- Wake geometries, as are important for 
analysis of wake effects on blade loads, are not necessarily well 
displayed by plots of averaged velocity. If it is accepted that 
the present model is physically realistic (based on the average 
wake-induced flow comparisons) , then the distorted wake geometry 
should also be realistic. The existing measurements are not suf- 
ficient to allow "final" evaluation of the model's acceptability. 
Instantaneous velocity measurements would provide a much more 
complete data set for evaluation of wake geometry programs. Such 
measurements would be extremely difficult to obtain, of course, 
and do not exist at this time. Distorted wake geometries for the 
three advance ratios for which flow calculations were made are 
shown in Figures 24 through 26. In the x-z plots, (the side 
view, looking from port to starboard) the z-coordinate scale is 
expanded 5 times. The plots show only one tip vortex, which is 
drawn as a continuous curved element. The computer model used 
straight elements and has twelve azimuthal steps per revolution 
(for these cases). As is seen from these plots, there is con- 
siderable wake distortion, much more in the lower advance ratio 
than in the higher advance ratio cases. The importance of the 
wake distortion depends on what is being studied, flight condi- 
tions, rotor configurations, etc. For example, the plots of 
Figures 24 and 25 indicate that the rotor is operating such that 
it may pass through its own wake, for these particular flight con- 
ditions and this rotor configuration. There is not any blade- 
tip vortex intersection for this particular wake-rotor position 
combination, but it appears that such intersections do occur where 
the wake has been distorted up through the rotor disc. Wake dis- 
tortion effects could be important, it appears, for many rotor 
combinations and flight conditions. Fixed skewed helical wake 
geometries would not allow prediction of such blade-vortex inter- 
sections as the present model appears to predict for single rotor 
helicopters, for example. It should be noted that the wake geom- 
etries shown in Figures 24 through 26 are for lightly loaded 
rotors, with values of C T on the order of .0035 and for disc 

loadings of approximately 48 newton/meter 2 . Considerably higher 
disc loadings and thrust coefficients are common on helicopters, 
and may result in higher downwash with correspondingly less chance 
for blade-wake intersections. 

Instantaneous velocity measurement data is not known to be 
available, so no comparisons on that basis may be made at the 
present time. However, instantaneous velocities at two positions, 
as a function of time, were given in reference 10. Figure 27 shows 
a comparison of the measured average velocity, computed average 
velocities, and computed instantaneous velocities using the single 
vortex model of reference 10 and the multiple vortex model dis- 
cussed herein. The calculated v(t)'s are more similar for z = 

-0.07 than for z = - 0.2. The v(t) calculated at z = - 0.07 
with a single tip vortex shows a slightly smaller range of values 
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than the v(t) calculated with the system discussed herein; the 
phasing of the largest variation is practically identical, and 
the overall shapes are similar. The v(t) calculated at z = - 0 . 2 
with a single tip vortex shows a much wider range of values than 
the v (t) calculated with the more realistic wake model, which 
indicates that a vortex passed close by z = - 0 . 2 in the tip- 
vortex-only model, but was not as close in the full mesh- 
modified wake model. The single-tip vortex model gives aver- 
aged velocities on the order of 20% higher than the measured 
average velocity, while the model presented here gives an aver- 
aged velocity less than 1% lower than the measured average 
velocity. As may be seen from Figure 19, there are a few 
points at which the single tip vortex predicted v more closely 
then the present model, but values of v(t) for other points were 
not presented in reference 10, so comparisons at other points 
cannot be made. 

Vortex blade interactions .- The possibility of vortex-blade 
loading due to radial motion of the point of interaction was dis- 
cussed in reference 18. Such problems led to suggestions for 
helicopter rotor systems with shaftwise separation of rotors, 
rotors having different lengths (and other physical differences) , 
and nonuniform azimuthal blade spacing, in hopes that such systems 
may be useful in alleviating blade-vortex interaction problems. 
Wake geometry and wake-induced velocity influence coefficients 
were computed for several such rotor systems, and wake geometry 
plots for some of those systems are shown in Figures 28 through 
33. Only two of the four tip vortices are shown, to reduce the 
complexity of the plots . The vortices are shown as continuous 
curves, but are based on sixteen steps per revolution. All co- 
ordinates are nondimensionalized by the reference rotor radius. 

Table 6 shows the various model parameters used for the wake 
runs. Figure 28 shows the free wake for Case 2, which is a stan- 
dard four-bladed helicopter configuration. The model used for 
this case was one with two two-bladed rotors with no shaftwise 
separation, 90° azimuthal spacing, advance ratio of 0.2, and with 
rotors rotating in the same direction. Figures 29 and 30 show 
the free-wake geometry for two rotor systems each with shaftwise 
separation of one chord for Cases 3 and 4 , which are for phys- 
ically identical blades, advance ratio of 0.2, and where Case 3 
has an azimuthal separation of 90° and Case 4 a separation of 
30°, respectively. Figure 31 shows the geometry for Case 7, 
which is for a system the same as shown in Figure 30, with an 
advance ratio of 0.3. Figures 32 and 33 show geometries for 
rotor systems with rotors of different lengths. Figure 32 pre- 
sents the geometry for Case 10, which is for a system with an 
advance ratio of 0.2, a one chord shaftwise separation, a 30° 
azimuthal separation, and with the upper, forward rotor length 
being 0.85 of the reference rotor length. Figure 33 shows the 
geometry for Case 9 , which is for a system with an advance ratio 
of 0.2, no chordwise separation, 90° azimuthal spacing, one rotor 
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with blades having no airfoil inboard of 0.4 r, and the other 
rotor having a length of 0.7 of the longer rotor. In all cases, 
the rotor system cyclic and collective values were intended to 
result in a lift of 17,800 newtons and drag trim for a fuselage 
flat plate area of about 0.74 meters 2 . These values of collective 
and cyclic were input, as were values of blade coning and first 
harmonic flapping. For Figures 28 through 33 the shaft axis is 
perpendicular to the tip-path plane. 

There are both differences and similarities in the various 
wake geometries shown. The self-induced effects of the wake 
elements upon themselves due to, their gross geometrical arrange- 
ments and circulation distributions appear to be the dominant 
source of free wake distortions. These gross distortion patterns 
are similar for all advance ratios and rotor configurations. 
Generally, near the rotor disc the horizontal distortions appear 
to be less severe than the vertical distortions . The edge of 
the wake, near y = ± 1, tend to roll upward, and there are 
smaller upward distortions near y = 0 at both forward and aft areas 
of the wake. In the region either side of y = 0, vortices corre- 
sponding to the wake deposited near the leading edge of the rotor 
disc generally appear to pass from below to above the rotor blades 
as they move aft. The x-position where they move through the 
disc depends on advance ratio. One or more vortex-blade inter- 
actions could result from such distortions. The loops near y = ± 1 
generally extend from below to above the rotor disc and are inside 
the edges of the disc. These loops begin to form almost immedi- 
ately after the tip vortex is deposited. After approximately one 
revolution the side loops have begun a gross roll-up pattern, 
especially on the advancing rotor side. Except at low advance 
ratios, the portion of the wake having this roll-up pattern has 
been transported aft of the rotor disc, so is not a region of 
potential vortex-blade interaction for the rotor systems and 
flight conditions investigated herein. The free wake distortions 
become increasingly pronounced with decreasing advance ratio. 
Figures 30 and 31 are wake geometries for identical rotor systems 
with advance ratios of 0.2 and 0.3, respectively, and the computed 
distortions are much smaller for the 0.3 advance ratio than for 
the 0.2 advance ratio. 

There are some characteristics which vary from case to case, 
and which may influence the type or severity of vortex-blade inter- 
actions. The wake from two rotors of equal lengths with 90° blade 
spacing, as shown in Figures 28 and 29, indicate that the distor- 
tion of the tip vortices near y = 0 is approximately the same from 
each blade. That is, the vortices are approximately at the blade 
tip elevation near y = 0 approximately one revolution after they 
were deposited in the wake near x = - 1. In contrast to this, the 
wakes from rotors of equal length but with the upper rotor less 
than 90° ahead of the lower rotor show that this rotor configura- 
tion results in this region of the tip vortex from the upper rotor 
being transported lower than the tip vortex from the lower rotor 
in approximately the same time. Such distortions can be seen for 
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Cases 4 and 7 in Figures 30 and 31, respectively. The same types 
of distortions were also observed for Cases 5 and 6, where Arp was 
20° and 10° , respectively. This type of distortion is that which 
would be expected from the intertwining of two tip vortices , and 
is a direct result of the nonuniform azimuthal spacing. It is 
possible that such intertwining could result in different types 
or severity of vortex-blade interactions. (Such results were 
not observed from the subsequent blade loads runs, as will be dis- 
cussed later.) Wake distortions for systems having rotors of 
different lengths generally showed more severe distortions than 
for systems having rotors of the same length. Figure 32 shows the 
wake geometry for Case 10, which is similar to Case 4 (shown in 
Figure 30) except that the upper rotor is 0.85 as long as the 
lower rotor, and the chords on both rotors were increased to re- 
tain the same amount of blade area for the total system. The 
intertwining process appears to be present, but it does not appear 
to have progressed as rapidly as for the cases with rotors of 
equal length. This is probably due to the tip vortices being 
positioned further apart (shortly after they are deposited) , and 
thus resulting in less interaction. The wake geometry shown in 
Figure 33 is for Case 9. Tip vortex distortions are more severe 
for this case than for any other case with an advance ratio of 
0.2. This appears due to the interaction of the short blade's 
tip vortex with that from the preceding long blade, near y = 0. 

Two of the tip vortices have been omitted from Figure 33. The 
location of the tip vortex from the long blade which is not shown 
is such that interactions with the following short blade's tip 
vortex are possible. 

Calculations of wake-induced velocity influence coefficients 
for all cases were made with built-in limits on the magnitude of 
the off-diagonal elements of the influence coefficient matrix. 

The number of times in which the limit was used to reduce the 
size of the off-diagonal elements was small. The use of such a 
limit was found to be necessary in some cases during preliminary 
calculations in order to allow convergence of the blade loads 
iterative solution methods. Such limits were not felt to be over- 
ly restrictive, but could diminish the severity of vortex-blade 
interaction effects, and thus represent a limitation of the pres- 
ent model and calculated results. 

Wake geometries shown in Figures 28 through 33 correspond to 
only one of sixteen reference blade azimuthal positions. Wake 
geometry plots could also be developed for any or all of the 
other blade azimuthal positions, and the tip vortices for other 
blades could also be shown. It is felt that presentation of such 
plots would add little to the basic understanding of the wake 
geometry model obtained from the plots already presented. Figure 
34 presents the wake for all four blades for Case 4 (corresponding 
to Figure 30) . The tip vortex intertwining process is observed 
for all four vortices. The vortex elements from the retreating 
blades which have been deposited for approximately one-half 
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revolution have approximately the same elevation at the front 
near y = 0 , but with the originally "upper" vortex now slightly 
lower than the originally "lower" vortex. It appears that the 
intertwining process does not continue after the initial exchange 
of positions, and that the vortices near y = 0 which were depos- 
ited near x = - 1 (forward part of rotor) may in some cases be 
sufficiently distorted to cause subsequent vortex-blade inter- 
actions near x = + 1 at the rear of the rotor. Such interactions 
appear to be likely, as seen from Figure 34, between the upper 
blade and the lower blade's tip vortex, near x = + 0.7 or + 0.8. 

In general, it does not appear that the various rotor systems 
provide a significant reduction in the number of vortex-blade 
interactions. However, the spacing (in time or azimuthal angle) 
and possibly the severity of the interactions are influenced by 
the various rotor systems. The present study indicates that some 
types of previously suggested vortex-blade interactions may not 
occur in the manner which would be expected for rigid wake models. 
For example, the type of traveling-force vortex-blade interaction 
suggested in reference 18 may occur, but the results of this 
study indicate that such' traveling-force interactions may be 
limited to a smaller portion of the radius than rigid wake models 
would predict. This does not mean that the resulting blade re- 
sponse would necessarily be milder. If this type of interaction 
occurred over only a limited radius range, but occurred near the 
blade tip where the dynamic loads are high and at a rate such 
that blade response was high, a smaller radial range of interac- 
tion could be significant. It does appear that vortex-blade 
interactions may occur for more flight conditions and over more 
regions of the rotor disk than would be predicted by a rigid 
wake model. 

Comparison of computed and measured average wake flows 
(reference 17 ) indicates that the present model is acceptable 
for average wake flow calculations, and implies that it is accept- 
able for wake geometry and wake-induced velocity influence coef- 
ficient calculations. The model uses a limited number of trailing 
vortices and azimuthal steps, which is important in obtaining 
numerical results at reasonable cost. The basic wake model is 
very flexible. It may be used for any two-rotor system which has 
equal rotor speeds on each rotor. It provides a basis for exten- 
sions to non-steady-state flight conditions and steady-state 
maneuver flight conditions. Running times are strongly dependent 
on the total number of elements generated. Wake flow runs for 
two and three revolutions of wake with four trailing vortices per 
blade in the modified wake and twelve azimuthal steps per revolu- 
tion ran on an IBM 360/65 for about 9 and 18 minutes CPU time, 
respectively. Calculation of wake-induced velocity influence 
coefficients were made (after considerable program streamlining) 
on a CDC 6600, using overlay features to keep core storage re- 
quirements below 70,000 octal; all runs were made with four 
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trailing vortices in the modified wake and sixteen steps per 
revolution. Runs of one and one-half, two and one-half, and 
three and one-half revolutions required CPU times of approxi- 
mately 2.5, 11, and 28 minutes, respectively. Thus, for moder- 
ate to high advance ratios, the present program and extensions 
of it would provide a flexible tool for use in further studies 
of one of the most complex problem areas associated with rotary 
wing aircraft operations - that of determining the free-wake 
geometry and its effects on blade airloads. 


Blade Loads and Response Results 

Blade loads were computed for twelve rotor systems - flight 
condition combinations, with more than one set of data being used 
for some cases . Output from the wake geometry program was used 
in all cases to provide starting values for bound circulations 
and wake-induced velocity influence coefficients. Blade proper- 
ties and model parameters as previously discussed, together with 
blade natural frequencies and mode shapes provided the additional 
model properties used in the calculations. The convergence 
criterion for the program was based on the relative change in 
bound circulation values from one iteration to the next. This 
tends to insure convergence of blade airloads. It does not in- 
sure convergence of blade torsional or bending moments and shears 
in all cases. Runs for all but two cases of the twelve satisfied 
the convergence criterion that the bound circulation relative 
magnitude change be less than 0.05. Cases 1 and 9 did not con- 
verge, according to this measure of convergence. 

Case 1 .- In Case 1, the circulations' minimum relative magni- 
tude change was about 0.16, which was sufficient to result in rea- 
sinably good airload estimates, but which was not sufficient to 
result in converged generalized coordinates. The predicted air- 
loads for Case 1 and the measured airloads for the approximately 
corresponding flight case are shown in Figure 35. Gross charac- 
teristics appear to be reasonable, i.e. steady and low frequency 
components are approximately the same at both radial stations. 

The detailed agreement is significantly better near r = 0.55 than 
near r = 0.85. The differences at both r's may be due to the dif- 
ferent climb rates, which was zero in the calculations and 350 
ft/min in the measured data. The measured airload at r = 0.85 
appears to contain effects of blade-wake interactions near ip = 225° 
and 315°. These localized variations were not predicted, and could 
be due to omission of climb rate from the model. Computed and 
measured blade shears and moments did not agree for this case, 
although steady flap and lag motions were approximately correct. 
Measured and calculated blade moments did not agree, possibly be- 
cause the generalized coordinates had not converged. Finally, the 
response formulation used for these calculations had the lead-lag 
damper effects included in the generalized force rather than con- 
tributing to the damping term in the equation of motion. This 
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permitted lead-lag variations which were large, but which coupled 
only weakly with the air loads. Thus, the lead-lag modes were 
diverging while the circulations were converging. The program has 
subsequently been modified so that lead-lag damping can be in- 
cluded in the damping term, and so that convergence depends upon 
convergence of both circulations and motions. 

Case 9 .- The model used for Case 9 calculations was that with 
a standard radius rotor with no airfoil inboard of 0.4R> and a 
rotor of length 0.7R at the same shaftwise elevation, and with 90° 
blade azimuthal spacing. The rotors' chords were equal, and were 
adjusted to keep the total blade areas the same as in other calcu- 
lations. Neither airloads nor generalized coordinates converged 
for this case. Several combinations of blade cyclic and collective 
angles and blade torsional damping were used in an attempt to ob- 
tain a converged, trimmed-out calculation. Program limitations 
on damping magnitude, o , in equation (28) , (over-critical damping 

was not allowed) restricted these efforts somewhat. Output from 
both preliminary calculations and from the non-converged runs in- 
dicated that the longer rotor was carrying approximately two-thirds 
of the total load. Flapping and torsional oscillations on the 
long rotor tended to become more severe as the iteration procedure 
between airloads and response progressed. It appeared that the 
expected lift (and corresponding drag and moment) were too severe 
for the system. 

Cases 2 thru 8, 10, and 11.- Blade loads for Cases 2 through 
8, 10, and 11 are given in Figures 36 through 44, in terms of the 
magnitudes of the harmonics of the root shears and moments. Gen- 
erally, the shears and moments corresponding to blade flatwise 

bending, V and M , respectively, had a large steady component, 
y z 

a first harmonic on the order of 30% to 50% of the steady, and 

very small higher harmonics. The shears and moments corresponding 

to the chordwise bending, V and M , respectively, generally had 

z y 

a large steady component, relatively large first and second har- 
monic components, and, in most cases, a rather large seventh har- 
monic component, with other harmonic components being relatively 
small. In all cases, the third through sixth harmonic torsional 
moment components were higher in proportion to the steady than 
for the other moments. For cases with identical rotors with 
various rotor arrangements, the variations of azimuthal spacing 
resulted in less extreme variations in blade shears and moments 
than did variations in advance ratio, as may be seen from Cases 4 
through 8. Converged cases with rotors of different lengths were 
similar except that the second harmonic component of the second 
rotor's chordwise shear and moment was extremely high. 

Beyond these generally similar characteristics each case has 
distinguishing details. As a check, Case 2 is a two-bladed 
rotor model of a single four-bladed rotor system and this gave 
almost identical airloads, shears, and moments on both rotors. 

This agreement between rotors is a measure of the numerical and 
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model accuracy. As seen from Figure 36, the numerical and model 
accuracy appears to be excellent. This does not prove that the 
model is correct. It does indicate that sufficient rotor revolu- 
tions were used in the wake program calculations to yield per- 
iodic wake effects, and that the blade response program gives 
consistent results for this system of "identical rotors". Case 
3 and subsequent cases in which rotor systems were non-identical 
showed differences for the different rotors, as would be expected. 
In Cases 3 through 6, these differences are due to shaftwise 
rotor separation or nonuniform azimuthal spacing. All these 
cases had identical advance ratios, cyclic and collective pitch 
angles and shaft angles on both rotors, as shown in Table 6. In 
some instances the harmonics of rotor loads differed by as much 
as a factor of two for the same harmonic on the different rotors. 
Usually the two rotors' loads were quite similar. Calculated 
differences for shaftwise separation of one chord, as seen in 
Figures 35 and 36 were small, with the upper rotor having slightly 
higher first harmonic loads, and with other load harmonics being 
slightly different for the two rotors. The differences between 
harmonic magnitudes for the two rotors generally increased with 
decreasing azimuthal blade spacing, Aip , as may be seen from 
Figures 37 through 40. Figures 38, 41, and 42 show loads for 
Cases 4, 7, and 8, where only the advance ratios were different 
and had values of 0.2, 0.3, and 0.1, respectively. The higher 
advance ratio case had the largest differences in harmonic load 
magnitudes for the different rotors and these occurred in the 
first harmonic loads. The third and fourth torsional moment 
harmonics were lower for the high advance ratio than for the 
intermediate and lower advance ratio. This may be an indication 
of fewer vortex-blade interactions for the high advance ratio 
case, which would be consistent with the wake geometries gener- 
ated for these runs. Finally, the second harmonic chordwise 
shears and moments increased with advance ratio, as would be ex- 
pected from the airloads. Results for Cases 10 and 11 are shown 
in Figures 43 and 44, respectively. Case 10 had equal collec- 
tive on both rotors. The long rotor thrust is approximately 
1.5 times the short rotor thrust, and except for the second har- 
monic loads, the long rotor had higher loading harmonics than 
the short rotor. Case 11 had approximately equal thrust on 
both rotors. Again, most other load harmonics, except the 
second, were larger on the long rotor. The unusually high shears 
and moments for the second harmonic of the short rotor 1 s chord- 
wise shear and moment for both Cases 10 and 11 are a result of 
aerodynamic excitation of the first chordwise elastic mode. 

The natural frequency for this mode for the short rotor is 2.04 
times the rotor speed. All cases had significant second har- 
monic chordwise shears and moments, as would be expected for 
helicopters in forward flight. However, the natural frequency 
of the first chordwise mode of the longer rotor was 1.51 times 
the rotor speed, so it was not aerodynamically excited. These 
exceedingly large loads could be reduced significantly with 
stiffness changes near the blade root so as to remove the 
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chordwise natural frequency from the second harmonic range. 
(Similarly, mass changes on the outer portion of the blade could 
be used to change this frequency from near 2 ft.) 

Case 12 .- Case 12 was run to try to determine the effects of 
flexible blade motion on wake geometry, and thereby on the wake's 
contribution to blade loads. All other cases, for example, used 
a wake geometry that was generated by blades with initially spe- 
cified motions.. Case 12 calculations were done instead by using 
blade positions as determined from Case 2 blade loads in a subse- 
quent wake geometry run. The blade positions at which wake- 
induced effects were calculated were not changed. However, the 
elevation at which the trailing vortices were deposited correspond- 
ed to the flapping positions of the flexible blade. These were 
determined in the Case 2 blade loads run. In this manner, higher 
harmonic blade motion effects were approximated in the wake geome- 
try program's computations. (These higher harmonic variations, 
however, were relatively small.) The wake-induced velocity influ- 
ence coefficients from this wake geometry run were then used to 
compute blade loads and response. The time histories of the air- 
loads on the reference rotor for Cases 2 and 12 are shown in 
Figure 45. The general trends were very similar, but there were 
differences over the aft half of the rotor disc. These differences 
were due entirely to differences of wake-induced influence coeffi- 
cients. It is apparent, then, that rotor blade position did have 
an affect on airload time histories, and it appeared that the aft 
half of the rotor disc was most affected, and that the higher har- 
monics of airloads were influenced more than the lower harmonics. 
Such results indicated that where wake-rotor interactions are 
important, it may be necessary to use more realistic blade posi- 
tions in wake geometry calculations than are permitted by defini- 
tion of coning and first harmonic flapping motions as given by 
the present performance program. Blade coning (wake position) 
was computed to be lower by the blade loads program than by the 
performance (blade position) program. It was felt that this was 
the dominant source of the differences in the resulting wake 
geometry and airloads. The general characteristics of the wake 
geometries were very similar, except that the blade coning angle 
was much smaller. 

Study of the time histories of the radial distribution of 
shears and moments provided additional information which could be 
utilized in this study. Figures 46 and 47 show radial distribu- 
tions of the blade flatwise bending moments as a function of azi- 
muth which were computed for the rigid rotor of Cases 2 and 7. 

For both cases, the moment values peak at the root, with an out- 
board peak for some azimuthal positions which is near the out- 
board peak of the third flatwise normal mode. The two rotors 
had essentially identical time histories for Case 2, but were 
somewhat different for Case 7. There were no severe peaks in 
the outboard flatwise moment distributions for any cases. There 
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were also no clear indications of traveling-wave-type interactions 
The use of only three flatwise inodes may not be sufficient to 
clearly display such a phenomenon, however. More radial load 
points and more flatwise bending normal modes may be required be- 
fore such a phenomenon could be displayed by such time histories. 


CONCLUSIONS AND RECOMMENDATIONS 


The following conclusions may be stated as a result of this 
investigation: 

1. The wake geometry, wake flow, wake-induced velocity 
model is a satisfactory model for advance ratios of approximately 
0.1 to 0.3. 

2. The wake program allows considerable flexibility for con- 
sidering rotor systems of different designs. 

3. The basic wake formulation permits extension to other 
steady-state flight conditions in addition to level forward flight 
As a real-time analysis, it also permits extensions to non-steady 
flight conditions. 

4. Wake program running times and storage requirements, 

for advance ratios of more than 0.1, are moderate enough to allow 
program usage in design problems. 

5. While the blade loads program provided reasonable re- 
sults for converged cases , some limitations were found where con- 
vergence was difficult to obtain for certain rotor systems and 
flight conditions. 

6. The blade loads formulation permits the program to be 
extended to other steady-state flight conditions in addition to 
forward flight, i.e. maneuvers. 

7. The blade loads program requires small running times and 
moderate storage. This allows the program to be easily and 
economically applied for design problems. 

8. Calculated time-average wake flows were satisfactory. 

Recommendations with respect to two areas are suggested. 

The blade loads program as presently formulated has limitations 
which restrict its usefulness and generality. The following 
changes are suggested: 

1. Develop and program solutions to the governing equations 
for damping values which exceed critical damping. 
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2. Develop and program solution of the governing equations 
for the generalized coordinates as a coupled set of equations. 

As a minimum, damping coupling should be provided but preferably 
complete coupling. This means permitting use of non-orthogonal 
modes . 


3. Remove linear damping effects from the forcing function 
and include these effects as damping terms in the governing equa- 
tions . 

4. Introduce in-plane wake-induced velocity effects on 
blade airloads in both the blade loads and wake geometry pro- 
grams . 

Recommendations for further research are suggested as fol- 
lows : 


1. Generalize the wake program to simulate transient as 
well as steady-state flight conditions. 

2. Study the effects of a freely deformed wake on airloads 
and blade response under transient conditions including gusts 
and transition from one steady-flight condition to another. 

3. Investigate the simulation capabilities of the wake 
program at low advance ratios (to hover) , and establish the pro- 
gram's capabilities for such flight conditions. 

4. Investigate simulation capabilities of both the wake pro- 
gram and the blade loads program at higher advance ratios (above 
0.3) . 


5. Generalize the wake and blade loads programs to permit 
simulation of main rotor-tail rotor interaction effects and 
the helicopter fuselage .and/or tail fin surfaces. 
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APPENDIX 


Equations Programmed for Computation 


Equations which were programmed for computation, and which are 
a result of equations discussed in the analysis section of this 
report are given here for completeness. 


Equations Used in Wake Calculations 


Blade flapping motions. - Blade flapping motions are obtained 
from the solution of the equations for £ lf B 2 , and 3 3 in subroutine 
BETAS , where 


3 i 

' 1 

Ki - M 10 

" M n /2 

0 

1 

1 — * 

M 0 0 

e 2 


- Mu 

K 2 - m 10 

ftM 2 0 " Mj 3/2 


M 0 1 

^3 

j 


- m 12 

- ftM 2 0 “ Mj 3/2 

k 2 - m 2 0 


Mq 2 


(Al) 


where 


M 0 o = PiV 2 A 0 /2 + p 2 [aV(a 1 + a g ) - flw ± - A R V 2 /2] 
+ p 3 ft 2 a 0 - r 5 n 2 A R , 

M 0 i = p 1 (3V 2 a 1 /4 - Vw. + V 2 a ) + 2p 2 VA 0 

1 s 

+ p 3^ (fiai — 2A^V) r 

Mq 2 = PiV 2 a 2 /4 + p 3 ft 2 a 2 , 

M 11 = P 1 Vcft/2 , 

M 2 2 = — P 2^^ > 
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Mi 3 — — P i V 2 /2 


M 2 0 = “ P 3 ^ • 

K 2 “ Mj o = Kl ~ Mi o - ft 2 Io , 

K i “ M i o = k + u z (I b + m j_ > 2 Cj _ > 6 - p 2 c/ 2 ) 
P ± = r i+1 - for i = 1,2,3, 

r i = upcCR 1 - r o ) /i , 

A R = a R/ (R " r o) > 

A 0 = a 0 + a R r 0 / (R - r 0 ) , 


a R = total blade twist, radians, positive for decreasing 
angle of attack with radius. 


a ]_ / a 2 = amplitude of lateral and longitudinal cyclic pitch, 

respectively , radians , 

and all other quantities are as defined in the list of symbols. 

Wake-induced velocities .- The computational equations for the 
velocity components induced by the blade and vortex elements are 
given in reference 19, and are included here, with minor changes 
in notation, for completeness. All quantities used in the compu- 
tation of wake-incuded velocities are in nondimensional form. 
Lengths are nondimensionalized by the radius of the longest blade, 
R; velocities by ftR; time by 1/^ ; and circulations , r , by ftR 2 . 


The points in a nondimensionalized reference system of loca- 
tion of trailing vortices from the blades, given by x^, y^, and 

z ^ j are determined by the blade flapping positions, azimuthal 

blade position, and input control numbers. If the vortex element 
end point P.. has coordinates (x.., y.., z..) and if v (x,y,z,^), 

1.J XJ lj XJ X 

Vy (x,y , z , \p) , v z (x,y,z, i)i) denote the dimensionless fluid velocity 
components at the point (x,y,z) corresponding to the rotor at 
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azimuth position 41 , 


then wake displacements are given by 


X i+l,j ( * + A*) 

*1+1, j ( * + 

Z ±+lrj ( ' P + AIP) 
for i = 1,2,.. 


= x ± j (4>) + 

= y ± j (*> + 

= Z ± j (*) + 

NW and 


v x (x i j ’ 


V (x. . , 

y J-D 


v (x . . , 

Z 1] ' 


3 = 


y ij' 


y ij' 


i j 

1 , 2 , 


ZjLj 

Z £j f'/') Al|) 
z i j ,4>)A'I> 

. . ,NGJ 


(A2 ) 


where NW is the number of trailing wake elements per blade and NGJ 
is the number of trailing vortices from all blades; and Ai|> = 2ir/NA, 
where NA is the integer number of evenly spaced azimuthal stations 
into which the rotor plane is divided. The fluid velocities are 
given by 


v = V + p cos a 
xx s 


v = 

y 


V 

y 


v = V + y sin a 
z z s 


(A3) 


where V , V , and V are contributions from the wake, y is the ad- 
x y z 

vance ratio, and a is the shaft tilt angle. 

s 

The strengths of vortex elements once formed are assumed to 
remain unchanged as the wake developes, therefore only subscripting 
changes need be done, i.e.. 


r. . (* + A*) = r. . (i>) (A4) 

-■-I) 1 1 , 3 

for i = 2,3,..., Nw , and j = 1,2,..., NGJ. The circulations 
r i j are determined as discussed in the analysis part, where the 

r's in equation (12) correspond to the I^.'s. The core sizes 
in the full mesh wake are assigned according to 

a . ■ (<|> + A*) = [L. j(*)/L..(* + AiJ,) ] 1/2 a. .(^) (A5) 

lj llfj 1 J 1 1 f J 
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for i = 2,3,.../ NW and j = 1,2,..., NGJ, which conserves the volume 
of rotational fluid. The core sizes at the blade (i— 1) are 
assumed to be of the form 

j = AKI*(AT ) 1/2 

with the constant of proportionality, AKI, specified by input data 
Here Ar is the spanwise difference of the blade ' s . bound r s, and 
represents the vortex element's circulation. Trailing vortex 
lengths, L. ., are given by 


L. . 
iD 




+ - yij )2+ - z ij> 2 i 


2i V 2 


(A6) 


while shed vortex lengths, ^ , are given by 

i/j+i " X ij )2 + (Y i,j+i “ y ij )Z + (Z i,j+l " z ij )2] / 


D. . = [ (x. . 

’ - 


(A7) 


and shed vortex core radii are computed as in equation (A5) but 
with Dj. replacing L . . , etc. Let q , q , and q be defined 
as J X ij y ij z ij 



q x = v x G 

9 



x ij x 




q = v G 

y . . y 

f 


and 

q z = V Z G 

Z ij 2 

/ 

(A8) 


where v 

x 


v 

Y 


- (y - - * i+llj > - - ^MYij 

* (z - z ij> (x ij - x i 

* < x - x ij )( *ij - sr i+l ,j> - ^ 


ir-ia -l+i. j’ ' <x ■ x ij )(z ij 


•yi+i.j’ , 

z i+l , j 1 , 
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G = 


r . . (r . . + r . J _. .) 

ij 13 x+ 1,3 


2 ’ r ij r i +> ,j [(r ij + r i+i,j^ - V 1 


and r.. = [ (x - x ..) 2 + (y - y <.) 2 + (z - z..) 2 ] 1/ ' 2 


ID 


ID 


ID 


ID 


If point (x,y,z) does not lie on a vortex element then vortex- 
induced velocity components are given by 


NW NGJ 

V (x,y,z) = l l q 

i=l j=l ij 


(A9) 


and similarly for V and V . If the point P is a vortex element 
J y z c rs 

end point with coordinates (x , y , z ) then 

xS ITS xS 


V x ^ x rs '^rs ' z rs^ 


NW NGJ 

l l q x 

i=l j=l X ij 

j^s 


NW 

+ I 

i= l 


q 


X . 

is 


+ q e 


(A10) 


X 


i^r-i ,r 


and similarly for V and V , where q , q , and q account for self- 

V z s s s 

2 x y z 

induced effects. If r>l and P is not on the "edge" of a wake, 

rs 


q = m F 
s x 

x 


q = m F 
^s y 

Y * 


and 


q = m F 
^s z 

z 


(All) 
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(the S factor, called R in reference 10, was omitted in error from 
the e terms) 
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rs 


1/2 


rs 


^ x r-i,s x r+i,s^ + ^ y r-i,s y r+i,s^ 


f = 


2S - 4s* - L|. 1>SJ 


2S 


+ /4S 2 - 


J r-i , sj 


+ (z r-l,s - z r+i,s> 2)1/2 


/L for L 2 , - 6 2 + L 2 

r- 1 , s r- 1 , s rs rs 


/L . for L 2 . > 6 2 + L 2 

r-i ,s r-l , s rs rs 
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2S 


2S 


- /4S 2 - L 2 


rs 


+ / 4 s 2 - 


rs 


/L for L 2 - 6 2 + Ii 2 

rs rs rs r-i,s 


/L for L 2 > 6 2 + L 2 

rs rs rs r-i , i 


If r=l, inclusion of blade-induced effects is accounted for 
by defining 


q S (X 1S' y is' Z 1 S > * (X 2B’ y 23' Z 2B ) 

X X 


r =0 

2 s 


similarly for q , and 
s 

y 


%J X 1B' y is' z is> = (X 2a , y 2s , z 2s > 


r e 0 

2 S 




B 


^ y r . £ n p . 

2ttR s j=A ^ *1 


(A12) 


where p = 


1 + /l+6 j' 


1 + /1+afj/ 


1 + /1+jif 


1 + /I+cT^ 


and where 61 = r g Aiji/dj 
3 2 


r s - D lj) 


«1 = 


3b /4d . 
s' ID 


“2 = 3hg/4 (d ! j - D 1;j ) 
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where r is the radius to the point P 

s is 

is the distance from P 


b is the blade semi- 
s 


chord at that radius, and d 


ID 


is 


and the 


furthest removed point on element 


ID 


If P is on the edge of 
rs 3 


a wake, so that, for example P and P are on the same 

J. / S"" 1 J7 S 

blade's wake, but P is on another blade's wake, then 

3T / St 1 


>s x (x rs' y rs’ 


z ) = q (x , 

rs' s x r,s-i' 


■r,s-i ' 


Z r f s-i ^ 


r =o 

rs 


and similarly for q g and 


notation q (x,y,z) 
s 

x 



q . In equations (All) and (A12) the 
z 

means that this is the value 

0 


q (x,y,z) would have if r. .were set to zero. Thus, the curvature 
s i n 

X J 

is determined from the three nearest end points, and the induced 
velocity is determined by the adjoining vortex elements only. This 
type of calculation is necessary at both inside and outside edges 
of the wake, as well as at the blade. 


Modified wake parameters .- The method used to specify the 
locations, strengths, core radii, and velocity influence coeffi- 
cients is flexible and consistent, but is not the only possible 
method. Implementation of the method used is done within sub- 
routines, and different methods may be implemented by programming 
and using different subroutines. The positions of the modified 
wake element end points which "attach" the modified wake to the 
full mesh wake are determined by the matrix relationship 


[xm] = [A] [x] 

where xm is a column of NTVM x-coordinate elements, 

A is an NTVM by NTV matrix, and 
x is a column of NTV elements , and 
similarly for ym and zm. 

The A is the same for all blades. Ordinarily NTVM < NTV, and A 
has nonzero terms only on the diagonal and just off the diagonal, 
with one or two nonzero elements in each row. In this case the 
modified element end points are on a straight line between two of 
the elements of the full mesh wake. 

Circulation values are assigned in a similar manner, and give 
modified wake circulation values. I'm, of trailed vortex elements 
based on full mesh wake circulation values, r, according to 
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[rm] = [B] [r] 


where Tm is a column having NTVM elements , 

B is an NTVM by (NTV-1) coefficient array, and 
r is a column having (NTV-1) elements. 


The B matrix elements are defined by input. Modified wake ele- 
ment core sizes are defined by input, except for the tip vortices 
which are defined to be the same as the corresponding full mesh 
wake tip vortex core radii. 


Velocity influence coefficients for the full mesh wake, a, 
are conveniently used by the blade loads program, rather than 
velocity influence coefficients for a combination of full mesh 
and modified wake representations. The definition of a’s which 
is consistent with the definition of modified wake circulations 
is given by 

[ a ] = [B] T [am] 

where a is a column having (NTV-1) elements 

am is a column having NTVM elements, and 

T 

[B] is the transpose of [B] . 


Equations Used in Blade Loads and 
Response Calculations 

Blade and wake circulation equations .- The equations which 
were used to compute blade circulation are also in nondimensional 
form and may be written as 


r(r i , ^ j ) 

Ek 

i 




V ij ) 


1 / 2 . 


C 2 <0 ij> 




V. . c (a. .)]w. . 
13 l il 


U? . 
11 


V? . 

11 


^b ± { e, (r ± , ^ j ) + ie 0 (*j + 1 ) - e 0 (<l» j _ 1 ) ]/2 a*} 


- i 


U. .a (a . . ) + V. .c„ (a . . ) 
n 3-1 il Z 13 


U? . + V? . ) 
il il 


T7T 


l l b 0 (r . , ip . ) 
m n m mn 1 ' v i 


r (r ,41 ) 
m ,T n 


m 


(A13) 
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where 


w. . = M b o (r . ,\p . j 
13 i: i; m mn i' r V 

J m n J 


r Cr w ,.* ) 

m n 


m J 


(A14) 


the r's in the expression for w^j are from the previous iteration, 
c (n, .) is the lift coefficient, and is a function of Mach number, 

x» 

as well as angle of attack, a(a^j) is the slope of the lift coeffi- 


cient versus angle of 
number, 

attack 

curve , 

and is also a function 

of Mach 

U ij = r i + u cos a s 

sin ip j 

/ 


(A15 ) 

V i j = h(r ± , <bj) - y 

sin a 

s 

- w. . 

/ 

(A16) 

a ± j = © 0 j ) + <f> (r ± . 

*j> + 

0 l (r ± ) 

+ tan -1 (V. ./U. . ) , 

13 13 

(A17) 


and other quantities are as defined in the list of symbols. 

Airfoil characteristics.- The lift, drag, and moment coeffi- 
cients, and the slope of the lift curve are computed by analytical 
approximations to experimental data for the NACA 0012 airfoil, in 
the subroutine SERIES , as follows : 

c^ (- a,^) = - c^(a,M^) and is given for ranges in ot as 

0-a< . 22689 (1-M^) , c £ = 5 . 7296a//lHY^ 

. .29269 (1-M. T ) + (1.3M. -.59)a 

.22689 (l-M N )-a<. 34906 , c £ = ^ -£ — 

(.12217+. 22689M) /1^2 


. 34906-a<2 . 7402 , c„ = (.080373 sin a + 1.04308 sin 2a 

z J 1 

A-m2 


- .011059 sin 3a + .023127 sin 4a) 
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2. 7402 - a <3. 0020 , c % - (.4704 + .10313a) 

3 . 0020-o-tt , c £ = (- 17.550 + 5.5864a) /G? » CA18) 

c^(-a) = c d (a) and is given for ranges in a as 
0-a<. 22689 (1— M^) , c d = (.0060 + .13131a 2 ) 

.22689 (1-M^jj) ~a~v , c d = - 1 (1 . 1233- . 029 894 cos a 

- 1.00603 cos 2a + .003115 cos 3a 

- .091487 cos 4a) (A19) 

c m (-a) = ~ c m^ a ^ ’ eva l uate< ^ for (chord/2), and is given for 
ranges in a as 

0-a<. 22689 (1-M^) , c m = 1 . 4324a//l^J^ 

< .29269 (l-M^)+(1.3M^-.59)a 

.22689 (l-M N )-a<. 34906 , c^ = - y 

[.48868+. 907561^] 

. 34906-a<2.7402 , c = — (-.02827 sin a + .14022 sin 2a 

m I 1 

/l-M 2 

- .00622 sin 3a + .01012 sin 4a) 

2.7402- a<3.0020 , c m = (.4786 + .02578a) 
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< < 


3.0G20-a-TT , c = (-12.5109+3.9824 

m 


a) /J\~. 


% 


(A20) 


(-O/Mn) = aCafMj^) 


and for all ranges of a may be obtained from c^(a, r^) according to 

(A21) 


. 9 C 

a(c, 'V = 




The quantities are input from the wake program. The 

equation for the (r/b) quantities is solved by an iterative process, 
during which values of w^ ^ r etc. are recomputed as the 

sets of (r/b) quantities are recomputed. 

Blade loads quantities . - When the r 1 s have converged , blade 
load quantities are computed according to the equations 


F Z (r i'V _ u ii°ii ( “ii )u n + _i_ 


pfi 2 R 3 b . 


2 A ip 1 


b. 


r(r i'*i + l> _ r ^ r i^i-i) 


b i 


b i 


+ 2b. 


c (ct. ) Xj? . , 7 V? . ’ 

m 1,3+1 i,D+l 1,3+1 


- c (a. . ) Su? " + V? 

m' 1,3-1' 1,3-1 i,3“l 


+ c , (a . . ) V. . /u? . + V? 7 

d 1-3 ' 13 13 13 


(A22) 


F (r . , \l > . ) 
x l = 

p ft 2 R 3 fcK 


V. . c n (a . . ) u . . / 1 

- 1 1 ‘z 3 - - i a ♦ »d'«i j > o «' 6 n + v ii 

b i 


(A2 3) 


where u.j^= /u?T + V?7 
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Mp 

pfi 2 R 4 b? 


b : x 


2c Co . . ) (u? . +v 2 .) - srrlb. 
m 13 X] 13 oArf) } x 


. r(r i-Vi f 


b. 

x 


b. 

x 


3nb^ r 


A\p 


6 (r i / i|»j +1 ) - 26 Cr, 1|»_. 


3irb^ r 


2 A\l> 


Cr. ij>.) + 0 Cr . <K H 
tf 3 3 ~ l 

5(r i'*j + l > ' 5tr i , *j-l ) 


(A24) 


Blade response .- Conversion of the aerodynamic loads to a 
form suitable for response calculations is done by lumping the 
distributed loads at the mass points of the lumped parameter blade 
model. The distributed aerodynamic moment is transferred to 
quarter-chord, then the distributed loads are integrated, using a 
straight- line approximation between load points , to obtain lumped 

J_T_ 

loads at the mass points. The drag force at the i— mass point, 
for example, is 


f 


x 


i 


F' Cr) dr 

X 


(A25) 


R. 


where and R^ +i are midway between masses rru and nu + , and be- 
tween m^ +i and m^, respectively, and where Cr) is the straight- 

line approximation to the distributed drag load. Coordinate trans- 
formations are then applied which result in loads in the local 
blade coordinate system, and are given by 


F =-f cos - f sin 4> . , 

v i z i 1 x i 


F = f sin 4> . - f cos $. , 

w. z. 1 


M, = In + Z F 

d> . u . a. v. 
Y x x xx 


(A26 ) 


where z = distance of the elastic axis forward of quarter chord . 

cl < 

1 

The computation of the total force and moment at each mass, for 
every azimuthal and radial location is done according to 
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Q = mfi 2 sin 4> (h + e cos $) + mO le + h cos (4 - 0 0 ) / cos 4] 

V c 

-2mfiijie sin $ + , 

Q, = mfi 2 er - 2fi 0 [ (I + me 2 ) sin $ + inch. sine 0 /cos 4] 

ij) C X u 

- 2fi^ i (i x + me 2 ) + 2mfie(v sin 4 + w cos 4) 

Q w = mfi 2 cos 4 (h + e cos 4) - mh.0 c sin (4 - 0 o )/cos $ 

- 2mfieiii cos $ + F 

w 

Q, = - fi 2 sin 4 [ (I + me 2 ) cos 4. + meh] 

<f> XI 

- 0 [I + me 2 + meh cos (4 - 0 O )/ cos 4] 

C X 

+ sin $ (I + me 2 ) + M (A27) 

x © 


where the Q's, 4, <p , <j> , v, w and the F's and M are understood to be 
subscripted by i and j; m, h, e, r, and I are^understood to be 

subscripted by i, and 0 O and 0 and its time derivatives are under- 

c 

stood to be subscripted by j. Here, 0 c , the cyclic pitch is given 
by 


0 = A sin fit + B cos fit 

c c c 


h, the horizontal separation between the elastic axis and the pitch 

t h 

axis, say at the i— station, is given by 

i 


h 


i 



cos 


$ 


l- 


l 


l is the change in distance of the elastic axis forward of the 

pitch axis in the local blade coordinate system; 4. . is the total 

th 

angle between the rotor plane and the chord at the i — radial and 
j— azimuthal position, and is given by 


$ 


ij 


i 

0 o + I 

j A=1 
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where 0o is the nominal pitch angle, and is the increase in 

twist angle at the station. Also, the radius to the i— 

mass is given by 

i 


r . 

i 


5 


+ 


I £ 

n=i 


n 


t 


and other quantities are as defined in the list of symbols. The 
values of v, w, <f> , and ^ in the above equations are from the pre- 
vious iteration. To account for a lag damper, the following is 
added to Q and Q , respectively: 
v 2 W 2 

AQ = (k £ /Jl 2 )i 1 sin $ 2 
AQ W2 = (k A /£ 2 )^l cos $ 2 


where k^ is the lag damping coefficient, ft-lb/rad/sec. 

The generalized forces acting on each normal mode are computed 
for each azimuth, according to 


F (t) = y (Q A (K) + Q A (K) + Q , 1 
K i v i v i w i w i ^i 


(K) + (K) > 

1 


+ 1 

J 


ct kj c j 


(A28) 


where is from the previous iteration. 


, = _ c a (k) .(J) 

J KJ C D* % 


for J fi K 


KK 


2 ff K “K ' 
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o is the average damping coefficient for the K — mode 
(read as input) , and 



is the mode shape quantity for the "q" type of elastic 

x. X- ' 4 - 

deformation, at the i— radial location, for the K— — 
mode. 


The 0 KJ terms may be thought of as damping coupling terms . The 
damping term, o K , is presently defined in terms of the mode shape 
quantities as 


a 


K 




2 

/2» k 


Note that the average damping coefficient, a , occurs on both 

i\ 

sides of the governing equation for c„(t), so effectively cancels 

i\ 

with convergence. (It should also be noted that the present pro- 


gram is inoperable for a v > 1, as calculation of vl-a* then re- 
suits in a fatal error.) 

The governing equation for the K— - generalized coordinate, 
C K , is given by 


C K + 2o K u K + “k C K F K (t) 

The solution of this equation is obtained in integral form. The 
solution assumes periodicity of both forcing function and re- 
sponse. The generalized coordinate, z , and its first time 

derivative, z v , are given explicitly below. The second time 

derivative of the generalized coordinate, ? K , is obtained from 

the governing equation itself, after solution for and £ K< 


Then, and are computed from 

i rt -o to (t-x) 

C K (t) = — - I e A sin u> R (t-x) F R (t) dr 

J 
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(A29) 


"Vk* 


lo (s 2 +c 2 ) 
K l K K 


fVc K + V .-] 008 V 


C K I c S K^ S s ^ n w K t 

K * K K 


■ , ft (t-T) 

5 K (t) - - VkV« *J ' K 


COS “ K (t-T)F K (T)<}T 


'“kV 1 


<S K + °K ) l- 


- si + cl 

K C K K S K 


sin to t 

JA. 


+ C K X c 


■ s„I 

K K S K 


COS Q) t 


where 


U K W K ' * “K 


A- a l 


”°k“k T . 

s R = e sm io k T , 


“ a K a) K T 

c R = 1 - e cos Uj,T , 


fT -a to (T-t) 

I s “ I e sin u> k (T-t)F k (t) dx , 


ft -a to (T-x) 

x c ~ j e cos w K (T-x)F K (x)dx 

K 0 


T = 2 it/ Q , 


, (A30) 
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and the integrations are performed by an approximation known as 
Filon's rule, (See, for example, appendix of ref. 20.) Response 
variables are then computed from 


v(t> - z , 

♦ (t) = I A^ (K) (t) 

K 

and similarly for w, ip, y, and -J = 3v/3r) . 

Conversion of the response data to a form suitable for loads 
calculations is then done by computing loads input data by inter- 
polation of response_data. The mass points on either side of the 
aerodynamic points, r. , are determined, and response quantities at 
r^ are computed by linear interpolation from the mass points. 

Thus, dimensionless response quantities, i.e., the plunging velocity 

• 

h, the pitch 9, and slope £, are given by 
h = [v cos $ - w sin $ + z cos <|> + (z - h) (0 - 



+ 5v cos a T 

cos fit 

r 

(A31) 

+ 

o 

CD 

II 

CD 

6 + d> , and 

c v 



(A32 ) 

II 

I 

£ 

cos $ + ^ sin 

*) - 


(A33) 


These values of h, 0, and £ (computed at all aerodynamic radial and 
azimuthal positions) are used in the next iterative calculations of 
blade loads • 

Once convergence has been established between blade loads and 
response, blade shears and moments are computed in subroutine SHEAR 
according to 

T (t ) = l A t (K) S K (t) (A34) 

K 
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for the torsional moment, and similarly for the flapwise and edge- 
wise bending moments M z and and for the flapwise and edgewise 

shears V and V . 

y z 

Initial estimates .- During the first step of the various 
portions of the program which involve iterations due to non- 
linearities in the unknown quantities, estimates are used to im- 
prove convergence^ Subroutine WBAR computes an average wake- 
induced velocity w which is uniform over the rotor disc, and 
which is given by 

w = { . 5 (y cos a T ) 4 + c 2 - (y cos a T ) 2 ]} 1//2 (A35) 


where y is the advance ratio, c^ is the thrust coefficient and a T 

is the angle of attack of the rotor tip path plane. The first 
estimate for blade circulations for the blade loads program are 
input from the wake program, and are those values associated with 
the last steps of the wake program calculations. The subroutine 
START computes the first estimates for various blade motion quan- 
tities which are used in the first blade loads estimates. The 
plunging velocity, h, is a function of both radial and azimuthal 
coordinate, and is given by 


h = r(Kj sin \p - n i cos ip) 


(b/2)[(A c - < x ) cos \p - (B c + m) sin i|i] 


- y cos tj; ( k o + <i cos + n sin 


(A36 ) 


the torsional deflection 0 varies only azimuthally, and is given by 


0 = 0 0 + A sin ip + B cos \p , (A37) 

c c 

the slope, 5, also varies only azimuthally, and is given by 

5 = k o + K i cos ^ + ni sin , (A38) 

and the time derivative of the generalized coordinates, assuming that 
only the primary flapping mode (say k= 1) responds, are given by 
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CiO) = (fiR/A sin ip - m cosip ) 

N 

and (ip) - 0 for K>1 . (A39) 

Here A ^ is the tip flap deflection of the normalized mode shape 
V N 

for the primary flapping mode. In the above equations, tc 0 , k 1# and 
n j , are given by 

Hi = ri 2 / n 3 r 
where 

n 2 = (dp-l) (a 2 d 0 -diao ) ~ ( a ido“ c i a o) (d 2 /2-co) 
and 

n 3 = (d 0 -l) [d 0 (d 0 -l)- di/2] - (d 2 /2-c 0 ) [d 0 (d 2 /2+c 0 ) - d lCl /2] ; 

<l = (l/d 0 -l) [ai-Ciao/do - (d 2 /2+c 0 -d 1 c 1 /2 d 0 )m] 

and 

< o “ ( a o“^i n l /2 ) /dp • 

In these equations 

c o = kj/4 

c x = ykj/3 

d 0 = wg _ b l k 1 /3 

d l = - ubikj/2 

d 2 = v 2 k x /4 

a 0 = [e 0 (l+p 2 )/4 + yA c /3 - yB^/4 
- (w+ya T )/3 + 6 N (1/5+y 2 /6) ] kj. 
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*l = [(1/4+p 2 /8)B c + A c b x /3 ] k x 

a 2 = 1 2 vi 6 0 /3 + (1/ 4+3v 2 /8) A - B bj/3 

c c 

- yCw+ya T ) + 

bj = blade semichord at radial station 1, 
$ N = twist angle at tip, 

kj = 2TTpbR 4 /I 2 , 

Z & z = l + fili/ij, 

f R 

lx = j ra (r) (r-6 ) dr 
6 


I 2 


f R 

m (r) (r-6 ) 2 dr 

6 


m(r) = radial mass distribution per unit span, and 

are as defined in the list of symbols. The k T and o> 2 

i p 

to the blade loads program. 


other symbols 
are inputs 
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TABLE 1 


ROTOR BLADE MODEL PROPERTIES, CONTINUOUS FORM 


Blade property 

Numerical value 

Blade station range 

mass distribution 

6.2419 kilogram/meter 

0.0762 meters to tip 

torsional inertia 
distribution 

0.02593 kilogram/ 
meter 

0.0762 meters to tip 

flatwise stiffness 

1.9696 x 10 ^newton- 
meter 2 

0.0762 meters to tip 

chordwise stiffness 

111.977 x 10 4 newton- 
meter 2 

0.0762 meters to tip 

torsional stiffness 

3.9393 x 10 4 newton- 
meter 2 

0.0762 meters to 
0.5334 meters 


2.0816 x 10 4 newton- 
meter 2 

0.5334 meters to tip 

location of pitch 
and elastic axes 

center line of shaft 

0.0762 meters to 
0.5334 meters 

(coincident) 

0.085725 meters from 
airfoil leading edge 

0.5334 meters to tip 

location of center 
of gravity forward 

0.0 meters 

0.0762 meters to 
0.5334 meters 

of blade twist 

0.006858 meters 

0.5334 meters to tip 

blade twist 

none 

0.0762 meters to tip 
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TABLE 3. 



j Natural 
| frequency 
| w/fi , for 
ift=5 . 9167HZ 


1.081 

2.041 


3.006 


Tip deflection, 
per unit flap 
deflection 

w/v (j)/V 

m/m rad/m 


0.0611 -.074 

-19.9 -.115 


Mode type, as indicated 
by sketches 

] 1 — — ^ 

Flap,v i Edgewise ,w Torsion, <f> 



00 


Mode type name 

1st flapwise bending 
1st edgewise bending 

2nd flapwise bending 
3rd flapwise bending 

1st (control) torsion 
2nd edgewise bending 







•^1 





TABLE 5. 




















Figure 1 . 


Flow diagram of program useage 
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Figure 2 . Wake model with combination of "full mesh" 

wake and "modified" wake of trailing vortices 
only 
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First rotor is located at coordinate 
system origin. 




Second rotor location and 
rotation direction is determined 
by input. 

Figure 3 . Wake geometry and blade loads coordinate system 
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Chord = 1.125 ft 


Blade 1, standard blade, R = 210", r 0 = 21" 


2a) Chord = 1.125 ft 
2b) Chord = 1.227 ft 


Blade 2, tip radius = 0.85 standard tip radius, r 0 = 21" 


4 


Chord = 1.687 in. 


Blade 3, tip radius = 0.70 standard tip radius, r 0 = 21" 


Chord = 1.687 in 


Blade 4, airfoil root cutout to 0.40 tip radius, R 

r o 


210 " , 
84" 


Figure 8. Sketches of the four model blade 
configurations 



















Figure 13. Wake survey positions, as seen from plan view, 
in rotor tip path plane coordinate system 
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Figure 16. Comparison of computed and measured induced downwash 
below the rotor plane, p=0.095, C T =0.0032, x=0.5 
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Figure 17 . 


Comparison of computed and measured induced 
downwash below the rotor plane, u=0.14, 
C T =0. 00371, x=-0 . 5 
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V/V 0 V/if 


x = 0.5 
z =-Q07 


- - computed, "FULL" WAKE, RASA 

— computed, tip vortex only 
o measured, NACA DATA,NACATR 1319 


x = 0.5 
z *-0.2 


or' O 

oo 


Figure 19. Comparison of computed and measured induced 
downwash below the rotor plane, u=0.14, 
C T =0. 00371, x=0 . 5 



Figure 20. Comparison of computed and measured induced 
downwash below the rotor plane, u=0.232, 

C =0.00321, x=-0 . 5 














V /v 



Figure 23. Comparison of computed and measured induced 
downwash below the rotor plane, vj=0.14. 


C T =0. 00371, x=0 . 5 
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(c) rear view 

Figure 25. Single two-bladed rotor's tip vortex 

locations, p=0.141, a =-5.1° 

s 
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(a) plan view 



(b) side view 


Figure 26. Single two-bladed rotor's tip vortex 
location, m= 0.235, a =-8.5° 




















(a) plan view 



(b) side view 


Figure 33. Two two-bladed rotors' tip vortex locations 

p=0 . 2 , Az/c=0, A^=90° , a =-2.33°, R 2 /Ri=0.7 

s 


rotor 1 cutout to 0.4 Ri 








Case 2 

Az/c = 0 
Aip = 90 ° 
R 2 = Ri 
y = 0.2 



ade shear and moment magnitudes versus harmonic 
mber, m, for Case 2 




Case 3 

Az/c = 1.0 
A\p = 90 ° 
R 2 = Ri 



de shear and moment magnitudes versus harmonic 
iber, m, for Case 3 
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Figure 38. Blade shear and moment magnitudes versus harmonic 
number, m, for Case 4 











Figure 42. Blade shear and moment magnitudes versus harmonic 
number, m, for Case 8 


O o 




v 21.6 21.1 21.7 21.2 

20 , \ \/ / 


0 2 4 6 m 


T 

10 

0- 

Figure 44. 


m 


Case 11 

Az/c = 1.0 
Aip = 30° 

R 2 = 0.85 Rj 

P = 0.2 


Blade shear and moment mangitudes versus harmonic 
number, m, for Case 11 
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Flatwise bending moment, M , 

z 

radial distributions versus 
azimuth for Case 2, rotors 1 
and 2; u=0.2 
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radial distributions versus 
azimuth for Case 7, rotor 1; 

p=0 . 2 
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